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ABSTRACT 

This report describes the MITRE version of the TRACE-D nroaram 
now in operation. While the primary function of TRACE-D is orbit 
determination, options are also available in the program for trajec- 
tory prediction and observational data generation. A functional des- 
cription of these features is contained in this report along with a 
complete user*s manual and a brief program description. 



iii 



PREFACE 

The program described herein is a modified version of the TRACE-D 
program originally developed by the Aerospace Corporation in Los 
Angeles. Although the modifications made were many and included con- 
version of the program to FORTRAN IV and STRAP, restructuring of the 
program logic, and the addition of minor options and features, tne 
basic functions (i.e. , orbit determination, trajectory prediction, 
and observational data generation) remain intact. 

The people involved in various aspects of the effort, in addition 
to the uthor, were Dr. K. K. Maitra, G. M. Hyder, P. E. Shifres, and 
S. Schwartz. A great deal of appreciation is due to K. W. Hubbard, 
C. C. Tomes, and other members of the Flight Mechanic's Department 
of the Aerospace Corporation whose cooperation in obtaining the pro- 
gram and unlimited assistance throughout the task were invaluable. 

This report is intended solely as a user's guide to the TRACE-D 
program, and therefore contains only a brief functional description 
of the analytics. For the reader seeking the mathematical details 
and arguments the ideal reference is the original document^ issued 
by the Aerospace Corporation (see [1]). 



^ Material from this Aerospace document has been freely incorporated 
into this report wherever pertinent and is not to be taken as 
original . 
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Section 1. 
BASIC INFORMATION 



1. 1 Introduction 



This document is intended as a user's manual and therefore does 
not discuss the theory and mathematical details of TRACE-D except in 
a functional context. Section 1 of this report presents some program 
history, a brief summary of the modification process, and a statement 
of the program* s present status . Section 2 specifies the functional 
capabilities of TRACE-D while Sections 3 and 4 describe input and out- 
put (patterned after the original TRACE-D document) and Section 5 
contains programming information. Wherever pertinent, throughout the 
report, material from the Aerospace documentation has been freely in- 
corporated. 

1. 2 Program History 

The TRACE-D orbit determination program was developed by the 
Aerospace Corporation as a multipurpose, flexible computational tool 
for application to problems in orbit analysis. A major design objec- 
tive was that the program be able to provide answers to a wide range 
of problems relating to orbit and system design, space vehicle perfor- 
mance, and force model analysis. It was not and is_ not a "real time" 
program and therefore it is most effective for post-flight analysis 
and also for research and investigative purposes. The program has 
been used extensively and successfully by the staff at Aerospace in 
just such an environment. 

TRACE-D was originally written in FORTRAN II and FAP for use on 
the IBM 7094 computer in the "CHAIN" multi-coreload mode of operation, 
About thirty-percent of the program consisted of machine language 
(FAP) subroutines including the numerical integration package. The 
status of the program at the time of its release to us was reported 
to be operational for all functions. 



1. 3 Modification Process 

The primary objective of this effort was to produce a version of 
TRACE-D operable on the IBM 7030 computer here at MITRE. To accom- 
plish this, two separate yet simultaneous endeavors were carried out; 
one being to get the original TRACE-D program to run successlully on 
a nearby 7094 installation, and the other the actual conversion of 
routines to FORTRAN IV and STRAP and the reconstruction of the program. 
Upon completion of this phase, the next step was the checkout and 
"debugging'* of the 7030 version of the program. This task was greatly 
facilitated and accelerated by running test cases for all options with 
both versions (7094 and 7030) of TRACE-D and checking all results. 

I zing accomplished the primary objective satisfactorily, work 
began on minor adjustments and additions to the program. Such fea- 
tures as observation residual plotting, an expanded earth geopotential 
capability, and a dynamic constraint matrix capability for site loca- 
tion adjustments were implemented. At the present time the program 
is operational on the 7030 computer and has produced meaningful and 
significant results for many MITRE orbit-related problems. 



Section 2. 
MAJOR PROGRAM FUNCTIONS 

2 .1 Tralectory Prediction 

Basic to all the functions of the TRACE-D program is the genera- 
tion of a time history (either in a forward or backward direction; 
see Section 3) of the position and velocity of a space vehicle in an 
inertial frame of reference. However, this inertial ephemeris as 
well as its associated ground track and altitude history are, in many 
instances, of interest in themselves. Thus, these computations are 
performed and printed out in the trajectory prediction mode of TRACE-D. 

The motion of a space vehicle is generated by numerical integra- 
tion of the appropriate differential equations of motion. Using a 
Cowell formulation (i.e., expressing the total acceleration vector as 
three components in a cartesian system [x, y, z]) three non-linear 
second-order differential equations are derived, A predictor- 
corrector numerical integration procedure based on "eighth-order dif- 
ferences" is us< d to obtain position (x, y, z) and velocity (x, y, z) 
at any time t jased on the known position and velocity at some ini- 
tial time tQ. 

The total acceleration vector, mentioned above, is actually a 
sum of the effects of the perturbing forces that comprise the TRACE-D 
dynamic model. The earth geopotential effects (gravity) are included 
in the form of a spherical harmonic expansion with provision for 
zonal harmonics J2 through J15 and all tesseral and sectorial 
terms through J15 15. Effects due to other bodies in the solar 
system (Sun, Moon, Venus , Mars , Jupiter) are computed from inverse- 
square law formulas and positions of the other bodies are obtained 
from tabulated coordinates stored on magnetic tape (planetary tape). 
Acceleration due to atmospheric drag is assumed to be directly pro- 
portional to the square of the velocity of the vehicle relative to 
the air. The atmospheric density is obtained from one of two different 



model atmospheres that are incorporated into the program according to 
the user's choice. Instantaneous changes in the inertial velocity 
vector may be applied at specific times to simulate maneuvers such as 
orbit adjust or vehicle separation. Also, an included low-thrust ac- 
celeration term may be used to simulate thrust tailoff in cases in- 
volving large engines, long-term constant thrust, or, in some instances, 
leaking tanks or valves. 

Although all of the above effects are programmed into the TRACE-D 
equations of motion, the actual model that is employed may include 
some, all, or none of these perturbations. The user has the option 
to choose the model by setting a series of indicators. The compara- 
tive LSting of associated input and output quantities given in Table 
1 suggests the range of potential application. Complete instructions 
for preparation of trajectory mode input data and a sample of a typi- 
cal trajectory mode output listing are given in Sections 3 and A. 

2- 2 Observation Generation 

The object of the data generation function is to generate various 
forms of simulated measurements from a given definition of the ephem- 
eris of a vehicle and the location of observing stations on the earth's 
surface. These measurements may be any of twenty different types 
from up to one hundred different sites and may be generated at user- 
specified frequencies. Both "clean" and "noisy" data are available, 
the latter being generated by adding noise of specified variances and 
biases to the "clean" measurements. Upon user request, the visibility 
function (i.e., the rise and set times) of a vehicle with respect to 
one or more stations may be obtained without printout of any inter- 
mediate data. When data is generated, it may be written on magnetic 
tape or punched on cards and is formatted such that it is readily 
available as input to the orbit determination process. 

This option lends itself well to many applications, two of which 
are the following: studies of the effects of specific data imperfections 



on orbit determination convergence behavior and the fitting of data 
using a different force model than that used to generate the data 
giving insight into the real-world problem of fitting live data with 
programmed (and always less than perfect) dynamic models. 

Since the calculations performed for this option encompass all 
the calculations of the trajectory prediction mode, the input quanti- 
ties include all those listed in Table 1 plus the additional infoirma- 
tion in Table 2. Available output quantities are listed in Table 2 
also, but a more specific list of available data types with definitions 
is included in Sections 3 and A. 

2.3 Orbit Determination 

Orbit determination is the primary function of the TRACE-D pro- 
gram and encompasses both the trajectory prediction and observation 
generation modes. Stated in simplest form the orbit determination 
problem consists of extracting information from observations (possible 
20 types) of a space vehicle. The data is generally collected by a 
network of tracking stations on the surface of the earth. Nearly 
always, the information to be extracted includes the orbital elements 
but may include many other parameters of the dynamic model and of the 
tracking system. 

The process of extraction in TRACE-D takes the form of a general- 
ized least-squares differential correction procedure. More specifically, 
the goal is to determine values for a set of parameters (from the dy- 
namic and/or observational model) such that the differences between 
the actual input measurements and corresponding values computed from 
the model (usually termed observation residuals) will be minimized in 
a least-squares sense. The solution set may be comprised of not more 
than one hundred parameters, of which not more than sixty may be tra- 
jectory-related parameters (i.e. , initial conditions, geopotential 
coefficients, etc.). The amount of data that may be fitted is re- 
stricted only by the nature of the problem and not by the program. 
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Table 2. 
Input/Output Quantities Associated with 
TRACE-D Observation Generation Mode 



INPUT DATA 

1. All input listed in Table 1. 

2. Site location coordinates 

3. Visibility function indicator 

4. Observation types 

5. Observation frequencies 

6. Noise information 

7. Output mode (card, tape) 

8. Refraction specifications 

9. Uncertainty matrix for 

adBARV at epoch 



OUTPUT DATA 

1. Computed observations 

2 . Rise-set times (per 

station) 

*3. Observation uncertain- 
ties 



Is the propagation and effect of input 9 above. 



In fact, TRACE-D is capable of fitting observational data from six. 
vehicles in a single orbit determination procedure. It allows for 
estimation of the trajectory initial conditions and drag coefficients 
associated with five satellites in addition to the usual parameters 
in a single satellite run. 

The list of all possible input for the orbit determination mode 
is quite extensive as can be seen from Table 3. This is, of course, 
due to the comprehensive model programmed into TRACE-D for the pur- 
pose of wider and more varied applications. However, flexibility has 
not been sacrificed for completeness, so that the program is easily 
adaptable to either a non-real or more-real environment by simple in- 
put adjustment. The input and output capabilities of orbit determi- 
nation are listed in Table 3 and described in detail in Sections 3 
and 4. Table 4, which follows, is a complete list of all possible 
parameters that may be differentially corrected in TRACE-D, 
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Table 3. 

Input/Output Associated with 

Orbit Determination Mode 



INPUT 

1. All input in Table 1. 

2 . Parameters in solution set 

3. Initial conditions for 

satellites 2-6 

4. Bounds for parameter 

correction 

5. Stntlon locntlons 

6. Data errors and biases 

7 . Data editing criterion 

8. Refraction correction 

indicators 

9. Iteration count 

10. Constraint matrix for 

parameters 

11. Observations from all 

vehicles 

12. Tape specifications 



OUTPUT 

1. Chronological list of all 

observations by satellite 

2. Variational equation solu- 

tions 

3. Trajectory printouts 



4. Observation partial 

derivatives 

5 . Residuals 

6. Station and observation 

type RMS summary 

7. Computed and predicted RMS 

8. Corrections 

9. Covariance matrix 

10. Correlation matrix 

11. Solution set 

12. Tape of residuals 
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Table 4. 
Parameter List* 



1. Initial position and velocity components in either spherical, 
rectangular, or classical element form for up to six indepen- 
dent space vehicles. 

2. Reciprocal of ballistic coefficients for up to six different 
satellites. 

3. Up to six velocity increments (KICKS) for one satellite. 

4. Amplitude and time constant for an exponentially decaying 
low thrust. 

5. Zonal harmonic coefficients J2 through Jjq and all 
te ;eral harmonics ^2,1 through J6,6- 

6. Constant biases on all types of observations. 

7. Scale factors for range and range-rate observations. 

8. Time biases (i.e., biases in reported times at which 
observations were made). 

9. Latitudes, longitudes and altitudes above sea level of 
observing stations . 



On any one run up to 100 of listed quantities may comprise 
the solution set with not more than 60 from types 1-5. 
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Section 3. 
Usage 

3 • 1 Prel i minary Facts 

In order to run the TRACE-D pro.^ram the user need not concern 
himself with the internal elements of the program^ nor worry about 
handling a cumbersome deck of cards. The entire program is stored 
on tape and can be used simplv bv attaching an appropriate data deck 
to the small basic running deck supnlied by MITRE and submitting the 
job for execution. The basic running deck contains only tape usage 
information and the necessary instructions to trigger execution, 
therefore it is never necessary for the user to alter it in any way. 
Due to the varied tape requirements for each of the major functions, 
there are three versions of the basic running deck and the appropriate 
one must be used for each option exercised. 

The input to the TRACE-D program falls into two categories : 
one which follows the normal FORTRAN IV rules of input and the other 
which is interpreted by a snecial routine in the program. The former 
method is used for all observing-station data (only a fraction of 
total input) and needs no description, while the latter method i^ 
used for all remaining input and requires a brief explanation . 

The input routine accepts each piece of data from a data card 
in a field of eighteen columns in length with four distinct subfields 
as shown on the following page. 



^ The user who is curious or interested in the actual programming 
facts will find these concisely presented in Section 5. 
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FIELD ^ 1 

LENGTH 



COLUMN NO . — > 1 



VARIABLE NAME 
OR 
LOCATION 



10 



VALUE 



6 7 



16 17 



N^ 



N„ 
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1. Prefix Field - The prefix field is one column in length and can 
occupy columns 1, 19, 37, or 55. The mode of the input that follows 
is dictated by what is punched (or not punched) in the prefix field. 
The following are the permissible prefixes. 

Prefix Data-Modes 

blank Floating point number, number 

contains or implies a decimal 
point. 

I Fixed point number, no decimal 

point. 

D Hollerith infomnation, data is 

not a number but a sequence of 
characters which can extend into 
the exponent field. 

G Used in connection with integra- 

tion parar'Mers, basic constants 
and the errth model. It indicates 
that card contains a symbolic name 
followed by a comment. 

H Followed by 1 or 2 in column 2. 

This is simply a header card to be 
printed at various places on the 
output. The number 1 indicates 
columns 1-80 and 2 indicates 
columns 81-132 on the printout. 

2. Location Field - The location field is 5 columns in length and 
starts in columns 2, 20, 35 and 56. The information punched in the 
location field is either a symbolic name or a subscript. A symbolic 
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name defines a location in core where the data is stored. In the case 
of a subscript, it is used with respect to the immediately preceding 
symbolic name. Therefore, whenever a symbolic name is the name of a 
table, a subscript of 1 is attached to that name and those subscripts 
greater than 1 indicate storage relative to the start of the table. 
All names and subscripts should be left-adjusted within the field. 

3. Value Field - The value field, 10 columns in length, starts in 
columns 7, 25, 43, or 61. This field together with the exponent field 
contains the numerical or hollerith data which is to be input. In the 
event the data is an integer (prefix = I) the exponent field is ignored. 
If the data is a floating point number (prefix column is blank, not 
punched) the decimal point, whenever it is not punched within the value 
field, is assumed to be after the last digit. Whenever hollerith data 
is input (prefix = D) the exponent field becomes an extension of the 
value field thus allowing for a total of twelve characters in the 
hollerith quantity. All data in the value field must be left-adjusted. 

4. Exponent field - The exponent field, two columns in length, starts 
in columns 17, 35, 53, and 71. This field allows for a range of 
floating point numbers from 10"^^ to 10" • An exponent of - 10 or 
smaller must start with the letter J (i.e., J0 = - 10, J5 = - 15, 
etc.). As in the value field the characters in the exponent field must 
be left-adjusted. 

Although the input routine will accept four input fields per card, 
it is not necessary to use all four. One field used per card is just 
as correct since the rest of the card is ignored. 

In addition to the precise field format there are a few other 
rules of input preparation, 

1. The actual order of the data cards is almost immaterial, the only 
restriction being that all values with subscripts in the location field 
are placed in the input deck relative to the last previous symbolic 
name. 

2. Any input which does not apply to the immediate run or is to have 
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a value of zero need not be specified, 

3. In the case of two appearances of the same location sjrmbol, the 

last one to appear defines the effective input, 

A. The basic input units are feet, degrees, and seconds. Therefore, 

whenever the units are not stated then they are a combination of the 

basic set. 

With the preceding basic description and the details on input 
preparation that follow, the user should have sufficient knowledge to 
use the TRACE-D program successfully. 

3.2 Basic Data 

The .erm "basic data" is defined as the data that is common to 
all the principal TRACE-D functions. This includes the constants for 
the integration procedure, the necessary physical constants, function 
indicator, trajectory specifications (like epoch time, initial condi- 
tions) and force model inputs. 

3.2.1 Constants 

The constants constitute two lengthy lists of input and in- 
volve quite a bit of keypunching. Because of this, there is available 
a standard set, in card deck form, which may be used as is, or altered 
for specific entries. This facility reduces the user's card prepara- 
tion time and eliminates possibility of error. 

In the event that the constants must be punched they must be 
formatted as shown on the following page. 
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3.2,1.1 Lines 1-19 Standard Set of Integration Constants^ 



LINE # 


PREFIX 


NAME OR SUBSCRIPT 


VALUE 


EXPONENT 




Cols. 


Cols. 


Cols. 


Cols. 




1,19,37,55 


2,20,38,56 


7,25,43,61 


17,35,53,71 


1 


G 


INTEG 






2 


I 


1 


1 




3 


I 


2 


2 




2 4 


I 


4 


1 




5 


I 


5 


1 




6 




11 


1.00002516 




7 




23 


1. 


JO 


8 




26 


1. 




9 




30 


1. 




10 




31 


.015625 




11 




32 


64. 




12 




33 


1. 


- 7 


13 


I 


34 


4 




14 


I 


35 


1 




15 


I 


36 







16 




37 


.001 




17 




38 


0. 




18 




40 


2820.1763 




19 


I 


41 


1 





^ For definitions of the constants consult Table 5. 

^ Notice that there are entries missing in the INTEG list, like the 
3rd one, for instance. This is because these entries are either 
not necessary for all options or are used by the program and not 
available to the user. This is true in both constants lists. 
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3.2.1.2 Lines 20 - 46 Standard Set Physical Constants^ 



LINE # 


PREFIX NAME OR SUBSCRIPT 


VALUE 


EXPONENT 




Cols. 


Cols. 


Cols. 


Cols. 




1,19,37,55 


2,20,38,56 


7,25,43,61 


17,35,53,71 


20 


G 


C 






21 




1 


.43752691 


- 2 


22 




2 


.55303935 


- 2 


23 




6 


.5 




24 




8 


.43752691 


- 2 


25 




9 


1. 




26 




13 


3280.8399 




27 




14 


57.2957795 




28 




15 


20925738. 




29 




16 


332951.3 




30 




17 


.0122999 




31 




18 


.814979 




32 




19 


.107821 




33 




20 


317.887 




34 




21 


95.129 




35 




22 


23454.865 




36 




23 


3443.9336 




37 




24 


20925738. 




38 




25 


348762.3 




39 




26 


32.174 




40 




30 


348762.3 




41 




31 


1.5 




42 




32 


1.0471976 




43 




33 


3.14159265 




44 




34 


298.3 




45 




35 


300000. 




46 




36 


82505.922 





For definitions of constants consult Table 5. 
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Table 5. 
INTEG and C Constants Lists 

INTEG PARAMETERS OF THE TRAJECTORY INTEGRATION 



1 1 Formulation. 1 - Cowell (eqs. of motion), 2 - Encke 

(note: the Encke-T formulation is presently not 
available) .^.^ 

2 2 Differential equation subroutine. 1 - AMRK, 2 - COW 

3 angle in degrees for changing mean to true equinox 

_A 1 Sun (A-9) are selectors for other-body perturbations 

5 1 Moon If tape unit number at numb (18) is not zero, 

6 Venus the perturbations are included or omitted 

_7 Mars according as the selector is non-zero or zero 

Jupiter 

Saturn 

Sun (10-15) are interpolation scale factors for use 

using contents of C(22) , the nuir^ber of earth-radii 



8 





9 








C(22) 



11 ■ 


l.UUUU^bib 


Moon 




12 


C(22) 


Venus 




13 


C(22J 


Mars 





14 C(22) Jupiter 

15 C(22) Saturn 

_] 6^ rJTP-Jlar velocity scale factors. (not now used, and 

thru all zero) . 

21 

22 ( used internally - Encke) 

23 IE-10 E Bar, COW subroutine truncation error control para- 

meter. 

24 (used internally - Encke) 
_2_5 (used internally - Encke) 



26 


1. 


A, see COW writeup. 


27-29 




(not used) 


30 


1. 


Initial time step size. 


31 


.015625 


Minimum time step size (1/64) 


32 


64. 


Maximum time step size 


33 


lE-7 


Kepler equation convergence criterion. 



34 4 Ratio of Cowell step size to Runge-Kutta step size 

35 1 if 1, do not recompute perturbations for corrector. 

if 2, do. 

36 Flacr for perturbation computation 
(used if INTEG(35) = 1) 



37 .001 Least squares convergence criterion (relative) 

38 , Least squares convergence criterion (absolute) 

39 (not used) 

40 2820.1763 Speed of light (earth-radii/minute) 

41 1 If non-zero, round position and velocity vectors 
for acceleration computation. If zero, truncate 
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INTEG 



Table 5 (cont.) 
PARAMETERS OF THE TRAJECTORY INTEGRATION 



42 



Reserved for later use 



C CONSTANTS 

1 .0043752691 



28 

29 
30 



Earth rotation rate (rad/mln) 



2 .0055303935 Gm, earth gravitation constant (ER**3/inin**2) 

3 Option for speed of light correction 

4 B=l-E, relative semi-minor axis of ellipsoid 

(computed in CSET) 

5 B**2/A**2 = (1-E)**2 (computed in CSET) 

6 .5 Factor for decreasing bounds in L. S. solution 

7 2*E-E**2 (computed in CSET) 

8 .0043752691 Atmosphere rotation rate (rad/min) 

9 1. Earth radius 

10 Input N for N-sigma residual editor 

11 Input scale factor for N-sigma residual editor 

12 If input non-zero, go thru 1st Iteration twice 

for sums 



13 


3280.8399 


Feet/kilometer 


14 


57.2957795 


Angle conversion factor 


15 


20925738. 


A, earth radius in feet 


16 


332951.3 


Relative mass of sun 


17 


.0122999 


Moon 


18 


.814979 


Venus 


19 


.107821 


Mars 


20 


317.887 


Jupiter 


21 


95.129 


Saturn 


22 


23454.865 


Earth-radii/astronomical unit 


23 


3443.9336 


Nautical miles (6076.1155 ft) /earth radius 


24 


20925738. 


I/O distance conversion factor 



25 348762.3 I/O velocity conversion factor 

26 32.174 Go (used for CDA/W and thrust/W) 

27 Input parameter difference for trajectory 
differencing 

Input threshold for percent difference for 

trajectory differencing 
Constant for doppler rate 
348762.3 Ft/sec per Er/min - 



31 1,5 Factor for increasing bounds in L, S. solution 

32 1.0471976 Rad/min per deg/sec 

33 3.14159265 PI 

34 298.3 Reciprocal of E = ellipticity 

35 300000. Critical altitude (ft) 

36 82505.922 Peg/day per rad/min 



37-39 



Direction cosines of body axis for look angle 
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3.2.2 



Functions to be Performed 



LINE // 



47 


D 


ITIN 


1234 





Line 47 contains the ordered list of all functions to be 
performed during a single run. Selection of functions is governed 
in accordance with the following code numbers. 

Code numbers Functions 

12 ' > Orbit Determination 

3 > Trajectory Prediction 

4 > Data Generation 

Up to twelve functions may be selected in any one run, i.e., the 
program will perform twelve specified functions in one continuous 
machine job. The example above specifies that an orbit determination 
be carried out first followed by a trajectory prediction and then a 
data generation. 

3.2.3 Trajectory Specifications 
3.2.3.1 Lines 48-54 Epoch 













48 


I 


YEAR 


1963 




49 


I 


MNTH 


6 




50 


I 


DAY 


15 




51 




TZNE 







52 




HR 


12 




53 




MIN 


45 




54 




SEC 


15.5 





Please notice that the field labels are no longer specified 
(i.e., Prefix, Value, etc.) but they are, of course, implied 
by the divisions. 
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In the usual case, wherein the year, month, and day are input with the 
year positive, the X-axis is directed to the vernal equinoic. Alter- 
natively, if the year is input negative, the X-axis would be directed 
to the longitude of Greenwich. The hour, minute, and second entries 
refer to midnight of zone time. Greenwich Mean Time is Time Zone 0. 

3.2.3.2 Lines 55-61 Initial Conditions 



55 


I 


ICTYP 


2 




56 




IC 


126.1 




57 




2 


31.23 




58 




3 


89. 




59 




4 


14. 




60 




5 


22600114 




61 




6 


25117.3 





Line 55 indicates which of the ten IC types (1, 2, , 9, 0) are 

entered in Lines 56 through 61. The alternative ICTYP entries are 
characterized as follows: 

a. IC Type 1 Earth-centered inertial cartesian coordinates 

(x, y, z, X, y, z in units of feet and feet 
per second) (see Section 3.1.1.1). 

b. IC Type 2 Spherical coordinates (a, 6, g A, R, v) in 

units of degr«es, feet, and feet per second 
(see Section 3.1.1.2). In Line 14, negative 
r is interpreted as height above the earth's 
surface in feet. In Line 15, if v is negative, 
cifcular velocity is computed and used. 

c. IC Type 3 Orbital elements (a, e, i, Q, w, t) in units of 

feet, degrees, and minutes (see Section 3.1.1.3). 

d. IC Type 4 Same as Item b above, with longitude X re- 

placing right ascension a. 

e. IC Type 5 No IC's input. The last trajectory point of the 

immediately preceding case is used. 

f. IC Type 6 No IC's input. The corrected initial conditions 

from the last previous tracking run are used. 



:22 



•^! 



g, IC Type 8 



h. 
i, 



IC Type 9 



IC Type 



Either Type 1 or 2 above, but in units of 
earth radii, minutes, and radians. Type 
number is entered at CPRAM (see Section 
5.1.2.2). 

Same as Type 1 above, but in units of 
earth radii and earth radii per minute. 

No IC's input. For a tracking run, two 
R, A, E sets are used from the data to 
calculate a set of initial conditions 
(see Section 3.6,6). 



3.2.4 Force Model Inputs 

3.2.4.1 Lines 62-68 Drag and T Matrix 



62 




DRAG 


.01 




63 


I 


2 


1 




64 




4 






65 




4 






66 




6 


-8.6 




67 




7 


-55.5 




68 


I 


TMATX 


4 





Line 62 contains the drag parameter CdA/W in square feet per pound 
and Line 63 contains the atmosphere model snecif ication (ARDC 59 or 
Lockheed/ Jacchia) . The ARDC model will be used when DRAG (2) is 
and the Lockheed/ Jacchia model when DRAG (2) is 1. If line 63 is a 1, 
Lines 64 and 65 contain quantities used in the Lockheed/Jacchia model. 

Lines 64 and 65 contain dj and d2» respectively. 
These are values of certain constants in the Lockheed-Jacchia atmos- 
pheric density expressions. (Consult reference 1 for details). 

An entry in Line 68, TMATX, will cause the T matrix 
(the 3x3 matrix describing the dependence of drag force upon vehicle 
position) to be used in the variational equations in accordance with 
the following options: 
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a. 


TMATX = 





b. 


TMATX = 


1 


c. 


TMATX = 


2 


d. 


TMATX = 


3 


e. 


TMATX = 


4 



T matrix is not used. 

Input 9p/9h is used with no earth 
flattening. 

Input 8p/8h is used with earth 
flattening. 

8p/9h is calculated with no earth 
flattening. 

9p/8h is calculated with earth 
flattening. 

Also, with TMATX non-zero, Lines 66 and 67 should contain input values 
for 8p/8h. Lines 66 and 67 contain 8p/8h for altitudes between 76 
and 108 n. mi. and between 108 and 376, respectively. 

It sho* d be noted that input values for 9p/9h must be used with the 
ARDC 1959 model. If a variable Cp term is desired, the drag table 
option may be utilized. In this case, the drag parameter may be con- 
sidered to consist of the product of two terms, (Cj)A/W) x c', wherein 
CdA/W is a constant which can be differentially corrected by the use 
of the variational equation and C/. may be considered a function of 
Mach number below a certain altitude and as a function of altitude at 
points above that altitude. Alternatively, C^ may be considered a 
function of time. In either case, use of the drag table is necessary. 

When the drag table is not used, input of CdA/W into 
the DRAG location does not change usual TRACE operation. In this 
case, Cj^ is automatically set equal to 1. If the drag table is used, 
the following additional inputs are required: 



a. 



DRAG (10) 



DRAG(ll) 



If DRAG(IO) = 0, tables not used. 

If DRAG(IO) = 1, Mach and altitude 
tables used. 

If DRAG(IO) = 2, time table only used. 

Altitude above which altitude table is 
used and below which Mach table is used 
(needed if (DRAG (10) = 1). 



DRAG(12) = Used by interpolation routine. 
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DRAG (13) 
through (35) 



Altitude table, or time table if 
DRAG(IO) « 2 (hi, C'(hi), h2, C^(h2) 
. . . , h , C^(h ), 0, 0, or 
ti, Cf)(ti), t2, Cfj(t2), . . . 



C^(tn), 0, 0) 
DRAG (36) « Used by interpolation routine. 



^n> 



DRAG(37) 
through (59) 



Mach table (stored as noted in Item c. 
above) . 



3.2.4.2 Lines 69-75 Other Body Perturbations 



68 


I 


CTAPE 


7 





If perturbations due to other bodies in the solar system are to be 
included in the trajectory calculations, a planetary coordinate tape 
must be mounted and the logical-tape unit number 7 must be entered 
at CTAPE. 



69 


G 


C 






70 




4 


1 




71 




5 


1 




72 




6 







73 




7 







74 




8 







75 




9 








Lines 69-75 contain the C constants that indicate which 
bodies are being considered. The six bodies corresponding to C(4) 
through C(9) are the Sun, Moon, Venus, Mars, Jupiter, and Saturn. 
Entering a 1 in C(4) implies that the Sun is to be considered, while 
entering a in C(9) implies that Saturn is ignored. Therefore, in 
the example above only the Sun and Moon are used in the other-body 
perturbation computation. 
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3.2.4.3 Lines 76-79 Exponential Thrust 



76 




THRST 






11 




2 






78 




3 






79 




4 







If an exponential thrust is to be used, the quantities T| in units 



of force/mass = ft/sec^, T2 in units of min"^ , and t^ 



and tf in 



seconds from midnight of epoch date must be input at THRST, THRST(2), 
THRST(3), and THRST(4) locations, respectively, 
3 . - . 4 . 4 Lines 80-96 Gravity Perturbations 



80 


1 I IFLAG 




81 


I 


28 


i 



Lines 80 and 81 above contain the input quantity 
IFLAG (28) = 0, This is an indicator that must be set in order to 
clear storage for the earth model that is to be specified. 



82 


I 1 lOB 1 10 1 


83 


I 1 2 6 


84 


13 10 


85 


14 6 



Lines 82-85 define the number of terms to be used in the 
earth geopotential expansion. lOB(l) is the number of zonals and 
I0B(2) is the index of tesserals. Therefore, in the example above, 
10 zonals will be used and tesseral terms up to and including the 
6,6 term. Lines 84 and 85 specify I0B(3) and (4) which are the number 
of terms to be used in the V-matrix. The V-matrix is a measure of 
the dependence of the gravitational force upon vehicle position and is 
included as one of the many necessary calculations for the observation 
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partial derivatives (for differential correction) . 

If these two values are set to zero the V-matrix calcu- 
lations are ignored. For a derivation of the V-matrix and its impli- 
cations Appendix B of Reference 1 is recommended. 



lowing order: 



86 




OBJZ 






87 




2 


1,826 


-2 


88 




3 


2,485 


-6 













Lines 86-88 specify the zonal coefficients in the fol- 



OBJZ (1) 
(2) 
(3) 



not used 
J3 etc. 



89 




OBJT 






90 




2 


1.556 


-6 


91 




3 


2.26 


-6 


92 




4 


1.72 


-6 












93 




OBLT 






94 




2 


16.18 




95 




3 


35.62 




96 




4 


108.2 





Lines 89-96 specify the tesseral coefficients and the 
lambdas to be used, in the following order: 
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OBJT 


^nm 


OBLT 


"^nm 


(1) 


J21 


(1) 


^21 


(2) 


J22 


(2) 


T22 


(3) 


J31 


(3) 


^31 


(4) 


J32 


(4) 


^32 


(5) 


J33 


(5) 


^33 


(6) 


ji+i 


(6) 


Ji+1 



etc. 



etc. 



3.2.5 Miscellaneous Inputs 

3.2.5.1 Lines 97-100 Instantaneous Orbit Adjusts 



97 




PKICK 


89020.31 




98 




2 


200. 




99 




3 







100 




4 


62. 



































Instantaneous orbit adjusts are input at PKICK. Line 97 contains ti 
(time of first orbit adjust, OAj) in seconds from midnight of epoch; 
Line 98 contains K (magnitude of velocity change of OAj) in feet 
per second; Line 99 contains 9 (yaw angle for OAj) in degrees, and 
Line lOG contains 9 (pitch angle for OA ) in degrees. 

The TRACE-D program will accommodate up to six orbit 
adjusts. Additional cards may be added as necessary for OA2 
through OAg. 



28 



3.2.5.2 Lines 101-107 Extra Kicks 



101 


I 


NXK 


3 




102 




XKICK 


63456 




103 




2 


.2 




104 




3 


88721 




105 




4 


10. 




106 




5 


101018 




107 




6 


.02 





Up to fifty fixed orbit adjusts (i.e., instantaneous 
changes of the in-track velocity component) may be input at XKICK. 
It should be noted that these orbit adjusts are not parameters for 
differential correction, but are applied in the equations of motion 
only and are independent of the PKICK inputs. The number of extra 
kicks (<_ 50) is input at NXK, and the table of times and AV values 
is input beginning at XKICK. The format is time, AV, time, AV, 
etc., in units of seconds from midnight of epoch day and feet per 
second, respectively. 

3.2.5.3 Line 108 Print Code 



108 


D 


PRCDE 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 



The print-code entry consists of two BCD words accommodating six 
character positions each. Entry of an X at any of these twelve 
positions will initiate corresponding outputs in accordance with the 
following: 

a. (1) Trajectory (trajectory only option). 

b. (2) Residuals (tracking only option). 

c. (3) Partials (tracking only option). 

d. (4) A^A after each iteration (tracking only option). 

e. (5) Variational equations (trajectory only option). 
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f. (6) Orbital elements (trajectory only option). 

g. (7) Do not print station locations (tracking and 

data generation) . 

h. (8) Not used. 

i. (9) Special trajectory prints^ 

j . (10) Orbital elements at ascending nodes only 
(trajectory only option). 

k. (11) Not used. 

1. (12) Suppress all trajectory print except ascending 
nodes (trajectory only option)^ 

3.2.5.4 Line 109 Equinox Precession Corrections 



109 




DALFG 


.0002 





The rotational position of the earth with respect to the inertial 
system is characterized by the right ascension of Greenwich at mid- 
night on the day of epoch (ag) and is computed with respect to the 
mean equinox of epoch date. An additional factor for correcting to 
true equinox of epoch date optionally may be input at DALFG in units 
of degrees. 

3.2.5.5 Lines 110-113 Direction of Integration 

The TRACE-D program is capable of generating a time 
history of position and velocity in a forward or backward direction 
in any one run. 

The inputs which determine this direction are INTEG(30), 
INTEG(31), INTEG(32). If these inputs are negative then backward 
integration from epoch time is assumed by the program. If positive, 



Prints will occur at times of maximum and minimum altitude above 
the oblate earth, at times when the flight-path angle equals 90 
degrees and at special latitudes and longitudes if values are 
entered (see Section 3.3.6). 

^ If the option to write a binary trajectory tape (B7) has been 
selected, the writing of that tape is controlled by the PRTIM 
entries (see Section 3.3.1). 
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as in the standard case, forward integration from epoch is performed, 



110 


G 


INTEG 




111 


30 




-1. 


112 


31 




-.01 


113 






-10. 



The example shown above indicates to the program that 
backward integration must be performed starting at epoch with an 
initial time step of -1. minute (see Table 5 for the explanation of 
INTEG constants), a minimum time step of -.01 minutes and a maximum 
time step of -10. minutes. 

When performing any of the three basic options, these 
inputs may be set negative. In trajectory generation a time history 
back from epoch will be produced. In the observation generation mode, 
observations previous to epoch may be generated. In the orbit deter- 
mination mode, data previous to epoch may be fit. Generally this 
makes TRACE-D a more flexible program by making the position of epoch 
and the known initial conditions less rigid. 

It should be noted that both backward and forward inte- 
gration can not be performed in the same run but two runs would have 
to be made, each with a different direction of integration. 
3.2.5.6 Line 114 Primary Scratch Tape 



114 



PTAPE 



PTAPE is the number of the primary scratch tape that 
the program requires and 4 is its only permissible value. This card 
must be included in every data deck. 
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3 • 3 Traiectorv-Only Data 

In addition to the basic data previously described, the following 
data is relevant to the trajectory generation. 



3.3.1 



Lines 1-7 Print Time Vector Option 



1 


I 


PRTIM 






2 


I 


2 






3 




3 






Ati 4 




4 






ti 5 




5 






At2 6 




6 

L _ 






t2 7 


7 















The above sequence of print times is for outputs selected by PRCDE 
entries (Line 108, Section 3.2.5.3). As many as nine sets of print 
intervals may exist (Line 2). In the case of the i^" set, output is 
from tj^_i to t-j^ at intervals of At^j^, with all times in minutes 
from midnight of epoch date if PRTIM = 1 or from epoch if PRTIM = 0. 
Additional cards may be inserted if 3 ^ n _£ 9 . It should be noted 
that a normal print at epoch is automatic. 

3.3.2 Lines 8-15 Variational Equation Partial Derivatives 



8 


D 


CPRAM^ ^ . 


















9 


D 


DPRAM 






















10 


D 


3 






















11 


D 


5 






















12 


D 


7 






















13 


D 


9 






















14 


D 


KPRAM 






















15 


D 


3 


i 

-U 


■ 
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An X entered in any CPRAM, DPRAM, or KPRAM character position causes 

the corresponding variational equation to be solved. Printout of the 

partial derivatives will occur only if an X is entered at Character 

Position 5 in the PRCDE print code entry location. The ordering of 

entries in the CPRAM, DPRAM, and KPRAM character position boxes is as 

follows: 

a. CPRAM (Initial Condition Parameters) (Line 8) 

The first position specifies which one of three types 
of initial conditions is applicable, and succeeding 
positions indicate the particular parameters that are 
desired in each case. Ordering of CPRAM parameter 
entries for initial condition (IC) Types 1, 2, and 3 
as shown below. 
(IC 
Type) 



(either) 


1 


X 


y 


z 


X 


y 


z 


to 













(or) 


2 


a 


6 


6 


A 


r 


V 


to 













(or) 


3 


a 


e 


i 


Q 


(A) 


T 


to 













Ordering of CPRAM IC Parameter Entries 

b. DPRAM and KPRAM (Differential Equation Parameters) 
(Lines 9 through 15) 

The ordering of DPRAM and KPRAM parameter entries is 
shown below, wherein Ti and T2 are the exponential 
thrust parameters and a-j^ , K^ 6y^ , and Qp^ are the 
OA number (1, 2, ..., 6), AV magnitude, yaw angle, 
and pitch angle, respectively. 

Sixty differential equation parameters is the maximum 
number which may be selected for any one run. 
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DPRAM j DRAG 


pi J 2 J 3 J'* i^S 1^5 


J? jJs 


J9 


Jio 


J2I 


|3 1 J22 


J±: IJ32 J 33 


Jhi JJh2 iJn3 


Jit it 


J51 


J52 


J53 


JSh 


I H 

5 i J55 


1 III 

■^51 ^62 ^63 J5'. 1^55 1-^66 i^2\ 


\22 


^31 


^32 


S3 


7 1 K, i^2 ks 


K. 


'^51 


S2 iS3 


Sit 


Ss 


Si 


S2 


S3 


9 ! ^54 1^55 ^65 1 Ua 


T. |t2 : 











KPRAM 



K. 



a, I K. 



>'2 i P2 



^P. r^ ! ^5 I ^y^ 



K. 



Ps 



i^e i Kg ; 



y3 



P6 



Ordering of DPRAM and KPRAM Differential 
Equation Parameter Entries 

3.3.3 Lines 16-18 Traiectory Comparison Options 



16 j 


1 
1 


IDIFF 


1 

1 


17 


I 

1 


NTAPE 


13 


18 


1 
j 1 


DTAPE 


1 


1 


1 ■ i 
1 \ 1 



Tape units and case indicators required for the trajectory differen- 
cing function are as follows: 

a. IDIFF = A regular trajectory run is indicated. 

b. IDIFF = 1 The reference trajectory will be written 

on the logical tape specified by NTAPE. 
If no entry is input at NTAPE, Logical 
Tape 15 will be used. 

c. IDIFF = 2 The differences between the present and 

reference cases are computed and written 
on the logical tape specified by DTAPE. 
If no entry is input at DTAPE, Logical 
Tape 14 will be used. The difference tape 
specified by DTAPE is rewound at the 
beginning of the case. 

d. IDIFF = 3 Conditions are the same as when IDIFF = 2 

except that the tape specified by DTAPE 
IS not rewound. 
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e. IDIFF = 4 The tape specified by DTAPE is rewound at 

the beginning of the case and unloaded 
upon completion. 

f. IDIFF = 5 The tape specified by DTAPE is unloaded 

upon completion of the case. 

The significance of the foregoing options is that if a single-compari- 
son case is to be processed, IDIFF = 1 is used for the reference case 
and IDIFF = 4 for the perturbed case. If a series of perturbed cases 
are to be processed, IDIFF = 1 is used for the reference case, IDIFF 
= 2 for the first perturbed case, IDIFF = 3 for all intermediate cases, 
and IDIFF = 5 for the last perturbed case. 

When this option is used, a special running deck must be obtained 
which defines the necessary tapes specified above. 

3.3.4 Line 19 Revolution Number 



19 




REV 


7 





If an initial value other than zero is desired for the revolution 
number, it may be input at the REV location. This value must be re- 
initialized for each individual case. 

3.3.5 Line 20 Trajectory Tape Generation 



20 



TTAPE 



If TTAPE is non-zero, a binary trajectory tape will be generated on 
Logical Tape 19 in accordance with the following input options: 

a. TTAPE = Tape will not be generated. 

b. TTAPE = 1 Tape will be rewound before generating 

but not unloaded after completion. This 
entry should be used for the first case 
when more than one case is involved. 

c. TTAPE = 2 Tape will not be rewound before generating 

and not unloaded after completion. This 
entry should be used for all intermediate 
cases. 
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d. TTAPE = 3 Tape will be rewound before generating 

and unloaded after completion. This 
entry should be used when only one 
case is involved. 

e. TTAPE = 4 Tape will not be rewound before generating 

but will be unloaded after completion. 
This entry should be used for the last 
case. 

This tape is not suitable for the trajectory-comparison option of the 

trajectory-only function. 

For this option also, a special running deck containing the proper 

tape definitions must be used. 

3.3.6 Lines 21-28 Latitude and/or Longitude Prints 



! "1 21 


I LATPR 


3 




22 




2 


10. 




23 




3 


15. 




24 




4 


20. 




n2 25 


I 


L0NPR 


3 




26 




2 


200. 




27 




3 


100. 




28 




4 


180. 





Line 21 contains n^ (n^ <. 10), or the number of special lati- 
tudes at which trajectory prints are requested, and Lines 22 through 
24 contain the special latitudes. Additional cards may be added if 
4 ^ n. :^ 10. 

Line 25 contains n2 ("2 < 10) . or the number of special longi- 
tudes at which trajectory prints are requested, and Lines 26 through 
28 contain the special longitudes. Additional cards may be added if 
4 <_ n2 i. 10. 

Note that an X must be entered in Character Position 9 of the 
PRCDE entry if either of the foregoing options are selected. 
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This concludes the description of trajectory related data. For 
data deck samples and running deck setups see Section 3.6. 

3.4 Observation Generation Data 

The following data is pertinent only to the data generation op- 
tion and is used in conjunction with the basic data (described in 
Section 3.3). 

3.4.1 Lines 1-4 Special Output Option 



1 




I FLAG 






2 


I 


6 


1 





If IFLAG(6) = 0, all generated data are printed. If IFLAG(6) = 1, 
rise, maximum elevation, and set times only are printed and the Data- 
Generation Specification Load Sheet II is not necessary except for 
listing of a card carrying TR in Columns 1 and 2 (see Section 3.4.5 
and 3.4.6). 



3 


I 


14 







If IFLAG(14) = 0, data are generated in time sequence until the 
available core space (bucket) is full. This output is then separated 
and printed by station in the same sequence as that of the input sta- 
tion cards. Further data are then generated until the bucket is again 
full, and the sort/print cycle is repeated. If IFLAG(14) = 1, data 
are printed as they are generated (i.e., in time sequence). 



4 




I 


ETAPE 


6 





If ETAPE is non-zero, a BCD radar observation tape will be gener- 
ated on the logical tape unit entered at ETAPE. The tape format will 
be the same as that of the tracking input data, including station 
locations and TF and TR cards. 
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3.4.2 Lines 5-23 Vehicle Attitude Specifications 



y 5 




YAW 


1. 




P 6 




PITCH 


90. 




r 7 




ROLL 


180. 















Vehicle attitude may be specified by inputting yaw, pitch, and roll 
angles in degrees in the manner shown above. These entries normally 
are introduced in conjunction with aspect angle computations. 



■ 
8 


YAW 






9 


811 


1. 




10 


812 


120. 




11 


813 


.01 




12 


814 


.1 




13 


815 







14 


816 


120. 




'' 


817 


121. 




16 


818 


-.9 




17 


819 


-9.0 




18 


820 







19 


821 


1420. 




20 


822 


1420.5 




21 


823 


360. 




22 


824 







23 


825 








The time history of vehicle attitude maneuvers may be specified by 
means of a table entered at YAW(811) . The format of this table, 
which is used in connection with generation of radar aspect angles 
and may consist of up to three sets of five entries each, is itemized 
in Table 6. 
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The yaw-, pitch-, and roll-angle values that make up these input sets 
of angular rates are assumed to change at the rate given over the time 
interval defined by the start and stop times. Nominal orientation is 
zero yaw, pitch and roll which corresponds to the condition where the 
vehicle body axis is normal to the geocentric radius vector, the nose 
of the vehicle is in the in-track direction, and the top of the vehi- 
cle is in the direction of the extended radius vector. 

Vehicle attitude at time of epoch and for the case where the en- 
tries in Table 6 are all zero (i.e., when nothing is input) is assumed 
to be the attitude specified at YAW, PITCH, and ROLL. If gaps in the 
time entries of the table are present, the angles are held constant at 
the last computed values. 

Table 6. 
YAW (811) TABLE FORMAT 



Entry 


Description 


YAW(811) 

YAW (812) 

YAW (813) 
YAW (814) 
YAW (815) 


Start time in minutes from epoch for 
first set of angular rates. 

Stop time in minutes from epoch for 
first set of angular rates. 

Yaw rate in degrees per minute. 

Pitch rate in degrees per minute. 

Roll rate in degrees per minute 


YAW(816) 

YAW(817) 

YAW(818) 
etc. 


Start time in minutes from epoch for 
second set of angular rates. 

Stop time in minutes from epoch for 
second set of angular rates. 

Yaw rate in degrees per minute. 

etc. 



3.4.3 Lines 24-38 Data Error Inputs 



24 



B NOISE 377777777 77 
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If NOISE is non-zero (a positive octal number) , normally distributed 
random noise with standard deviation and mean value specified by input 
at RAPAR is added to the generated data. The entry at NOISE is used 
to start the random-number generator. 



25 


M 


RAPAR 


05,99 












































•26 


D 


01,01 


H 


U 






A 


bIi 


A 


S 






27 




04,01 


.057 
























28 


D 


01,02 


A 


L 






H 


B 


I 


A 


S 






29 




04,02 


200. 

































The RAPAR array is used to indicate bias values for measurements 
to be generated. Line 25 merely indicates to the program that the 
array is of maximum size 05,99, Lines 26 and 27 specify that for 
station HU and the azimuth measurements the bias value to be added is 
.057 degrees. Lines 28 and 29 indicate that for station AL and the 
height measurement computed from AL the bias value to be added is 200 
feet. For additional entries the same two cards are repeated with in- 
creased subscripts . 
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30 1 




SIGMA 


200 




31 


t 
1 


2 


.1 




32 




3 


.1 




33 




4 


50. 




34 


I jISIG 


1 1 


35 


I |2 

1 


2 


1 


36 


I l3 

1 


3 




il 


I J4 1 112 i 



Lines 30-37 contain the observation data weighting factors. For 
each SIGMA entry, a corresponding entry defining the sigma set and 
data type appears in ISIG Lines 34-37. The ISIG entries are of the 
form 1001 + K, where I is the observation set number and K is the 
data type. Ten sets, corresponding to I = 0, 1, 2, . . . , 9 , may be 
entered. This selected value of I is the same as the entry In 
Column 5 of the station location card. (See Section 3 ,4 ,5 .) The 
data type, K, must be one of those listed in Table 7. 

In the above example, range, azimuth and elevation data would 
contain noise with a sigma of 200 feet, .1 degree and .1 degree res- 
pectively, while height would contain noise with a sigma of 5 feet. 
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REFR 



If elevation is a chosen measurement to be generated, then REFR 
is the input quantity that allows tor elevation refraction correction. 
The computed elevation, E, is altered to account for refraction, 
using either 

E' = E t Ogi cot E 
if E j^ 0.1 radian, or 

1 Hsi - 10^ 80 



E' = E + 



1000 12+iOOOE 6+lOOOE 
if E < 0.1 radian and n ^^ f 0. 
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REFR contains the ^g^ term, wherein i « 0, 1, 2, ..., 9. The 
appropriate value of i is entered on the station card in Column 6 
(see Section 3.4.5). Nominally Pso = 3.12 ^ 10"*^. Rise, set, and 
maximum elevation values are determined from the geometric E which 
represents the elevation before refraction correction is applied. Ad- 
ditional cards may be inserted at this location if necessary. 

Table 7. 
Data Types for ISIG Entries 



Data Type (K) | 



Data Description 



Symbol 



1 
2 
3 
A 

5 
6 
7 

8 
10 
11 
12 
13 
14 
15 
17 
18 
19 
20 
21 



Range 

Azimuth 

Elevation 

Topocentric right 

ascension 
Topocentric declination 

Topocentric hour angle 

Geocentric right 
ascension 

Geocentric declination 

Argument of latitude 

Cross plane 

Height 

X 

y 

z 

Range difference 

Range difference 

Range rate 

R difference 

R difference 



R 
A 

E 

a 
T 

6^ 
HA 

ag 
u 

V 

h 

X 

y 

z 

p 
Q 

R 
P 

Q 
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3.4.4 Lines 39-45 Orbit Covariance Matrix 



39 




ATAS 


1. 


-4 


40 




3 


.5 


-7 


41 




6 


2. 


-2 
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10 


4 


-2 


43 




15 


2000. 




44 




21 


1. 





If observation uncertainties are to be calculated, a covariance 
matrix for the ADBARV elements at epoch must be input in lower tri- 
angular form at ATAS. The order of the elements is Row 1/Column 1 at 
ATA£:(2), Row 2/Column 2 at ATAS(3). 

In addition, the following input is needed. 



45 


D 


CPRAM 


2XXXXXX 





The CPRAM variable indicates to the program that the covariance 
matirix has been input and in the ADBARV notation. (For a more com- 
plete explanation of the CPRAM variable see Section 3.3.2). 

3 . 4 .. 5 Station Cards 

The locations of tracking stations which are associated with 
the observations to be generated must be input by means of station 
cards carrying appropriate descriptive information. 

Figure 1 on the next page shows the proper fields for the in- 
put that is described below: 

a. Station ID Columns 1 and 2 (ST). 

No two stations may be identified by the same symbol. The 
three codes that are not permissible as station ID's are 
TS, TR, TF. 

b. Sigma Index Column 5 

The sigma index identifies the set of observation sigmas 
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.. ^ 



to be applied to data that is generated from this 
station (see Section 3.4.3). 

c. Type of Refractivity Correction Column 6 

This indicator selects the type of refractivity 
correction which is to be applied to elevation 
readings. (This input depends on previous in- 
put. See Section 3.4.3). 

d . Latitude Columns 9-17 

The latitude of the station is entered in 
degrees. 

er Longitude Columns 19-27 

The east longitude of the station is entered 
in degrees also. 

f. Altitude Columns 29-36 

The altitude of the station is entered in 
feet. 

g. P, Q Indicators Columns 38-39 and 41-42 

If P, Q, P, Q data (i.e., range difference and 
range rate difference between a master station 
and two slave stations) is to be generated these 
columns must contain the station ID's of the 
slave stations. Also, it is necessary to enter 
a station card for each slave station but col- 
umns 38-39 and 41-42 are to be left blank. 

The station cards must be preceded by two cards with END In 

columns 1-3, and must be followed by a card with TS in columns 1-2. 

3.4.6 Data Limit Specification Cards 

a. Station ID (ST) Columns 1-2 

b. Data Rate Columns 9-16 

Time interval in seconds at which data for a given 
station are to be generated and testing interval 
for rise set option. 

c. Minimum Elevation Columns 18-23 

The minimum elevation at which the vehicle is 
visible to this station. 



44 



d. Maximum Elevation Columns 25-30 

The maximum elevation at which the vehicle is visible 
to this station. Zero value implies a 90 degree limit. 

e. Maximum Range Columns 32-AO 

The maximum range in nautical miles to which the 
vehicle is visible. Zero value causes this test 
to be ignored. 

f. Start Time Columns 51-58 

Start time from midnight of start date in days, 
hours, and minutes. Zero value implies that 
epoch is start time. 

g. Stop Time Columns 60-67 

Stop time from midnight of start date in days, 
hours, and minutes. 

The set of cards must follow the TS card in the deck and them- 
selves be followed by a card with TR in the first two columns. 

3. A. 7 Data Type Cards 

The following set of cards (which are the last in the card 
deck) specify the exact type of measurements to be generated from each 
station. 

Figure 3 is a card image. 

a. Station ID Columns 1-2 

b. Observation Types Columns 7-33 

An X in the appropriate column will produce specific measure- 
ments according to Table 8. 
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NOTES 

ST Is the station Identification code. 

Station codes must be entered In Coluwns 38-39 and 
4i-A21f P, Q or P, Q data are Input. 

An asterisk {*) following a unit value Indicates 
that a decimal point Is necessary for that 
quantity. 
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NOTTS 

ST Is the station 
Identification code. 

Start and stop times 
are from midnight of 
epoch. 

If etart time Is 
epoch Itself Is used. 

An asterisk (*) fol- 
lowing a unit value 
Indicates that a 
decimal point Is 
necessary for that 
quantity. 
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NOTTS 

ST Is the Identifica- 
tion code. 

If P. P , Q, Q are 
selected the speci- 
fication of an ad- 
ditional station- 
Identification code 
Is necessary. 
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Table 8. 
Output Quantities Corresponding to 
Columns 7 through 33 of Data-Generation Type Card 



Column 


Output Quantity 


Unit 


7* 


Range 


n mi 

(ft on ETAPE) 


8* 


Azimuth 


deg 


9* 


Elevation 


deg 


10* 


Range rate 


ft/sec 


11-14* 


P, Q, P, Q 


ft/sec, ft 


15 


Azimuth rate 


deg/min 


16 


Elevation rate 


deg/min 


17 


Range acceleration 


ft/sec2 


18 


Mutual visibility (Output will be 
a list of numbers of stations 
visible at output time. Stations 
numbered in order of input on 
station cards. Number of stations, 
8 maximum) . 




19 


Geodetic latitude of vehicle 


deg 


20 


Longitude of vehicle 


deg 


21 


Surface range from station to 
subvehicle point 


n mi 


22* 


Altitude of vehicle 


n mi 

(ft on ETAPE) 


23 


Doppler rate 




2A 


Look angle (Angle between a vehicle 
axis and the station/vehicle line 
of sight. The direction cosines of 
the vehicle axis in the basic iner- 
tial system must be entered in C(37), 
C(38), and C(39). These quantities 


deg 



* These quantities are output on ETAPE 
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Table 8 (cont.) 



Column 



Output Quantity 



Unit 



25 



26 



27 



28 



29* 



may be input as constant or the user 
may provide a subroutine (FANG) to 
compute the direction cosines at each 
output point) . 

Observation uncertainties. (If in- 
verse A^A matrix for initial condi- 
tions is input at ATAS, the ADBARV 
elements are selected as parameters, 
and if an X is entered in Column 25, 
the [A^A]~^ is updated to observation 

times and the standard deviations in 

• • • 

the quantities R, A, E, R, A, E are 
derived and printed. The uncertain- 
ties are only those due to the un- 
certainty in the ephemeris which is 
implied by the given [A^A]~^ for the 
epoch conditions) . 

Angle kappa (K) . (Angle between 
station line-of-sight and geocentric 
radius vectors). 

Aspect angles. (Angle 1 ($) is de- 
fined as the angle between the vehicle 
yaw axis and projection of the station 
line-of-sight vector in the roll plane. 
Angle 2 (0) is defined as the angle be- 
tween the vehicle roll axis and the 
line-of-sight vector to the station. 

Signal attenuation = -40 logioR> where 
R is slant range in feet 

X, y, z (Same rectangular earth-fixed 
(X through Greenwich) geocentric quan- 
tities accepted as Type 5 observations 
for orbit determination) . 



Same units 
as obser- 
vations 



deg 



deg 



db 

n mi 

(ft on ETAPE) 



* These quantities are output on ETAPE 
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Table 8 (cont.) 



Column 


Output Quantity 


Unit 


30* 


Topocentric right ascension and 
declination 


deg 


31* 


Geocentric right ascension and 
declination 


deg 


32* 


Topocentric hour angle 


deg 


33* 


Vehicle-centered argument of 
latitude and cross-plan angle 


deg 



* These quantities are output on ETAPE 
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3.5 Orbit Determination Data 

This section gives a description of all the pertinent data for 
a differential correction run to be used In conjunction with the pre- 
viously described basic data. 

3.5.1 Lines 1-14 Multiple Satellite Option 

As was mentioned In Section 2, TRACE-D will accept measurements 
for up to six satellites on any one run. The method of specifying the 
input for these vehicles is shown below. 



Year 1 


I 


SAT2 


1964 




Month 2 


I 


2 


2 




Day 3 


I 


3 


10 




Hour 4 




4 


3. 




Minute 5 




5 


30. 




Second 6 




6 


52. 




ICTYP 7 


I 


7 


2 




IC 8 




8 


352. 




2 9 




9 


10. 




3 10 




10 


90.05 




4 11 




11 


165. 




5 12 




12 


22580632. 




6 13 




13 


25205.3 




DRAG 14 




14 


.015 





If observations for a second satellite are to be input, Lines 1 through 
14 are used for entry of epoch, initial conditions, and drag coeffi- 
cient for Satellite 2. Lines 1 through 3 contain the year, tnonth, and 
day, and Lines 4 through 6 contain the hour, minutes, and seconds, 
Greenwich time. Line 7 indicates the type of initial conditions that 
may be entered in Lines 8 through 13 (Type 1, 2, 3, 4, 8, or 9). Line 
14 contains the drag coefficient (CdA/W) . The input initial conditions 
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(IC) for satellites 3 to 6 are the same as those noted above except 
that the symbol SAT2 is replaced by SAT3, SAT4, SATS, and SAT6, as 
appropriate. 

3.5.2 Lines 15-31 Parameter Specification 

The method of selecting the parameter set for the differential 
correction process is described below. 

3.5.2.1 Differential Equation Parameters 



15 


D 


CPRAM 2 


X 


X 


X 


X 


X 


X 














u- 


16 


D 


CPRAM X 




























17 


D 


3 




























18 


D 


5 




























19 


D 


7 




























20 


D 


9 




























21 


D 


PSAT2 




























22 


D 


PSAT3 




























23 


D 


PSAT4 




























24 


D 


PSAT5 




























25 


D 


PSAT6 




























26 


D 


KPRAM 




























27 




D 


3 





























An X entered in any CPRAM, DPRAM, PSAT, or KPRAM character position 
causes the corresponding parameter to be differentially corrected. 
The ordering of entries in the CPRAM, DPRAM, PSAT, and KPRAM charac- 
ter position boxes is as follows: 

a. CPRAM (Initial Condition Parameters for Satellite 1) 
(Line 15) 

The first position specifies which one of these types 
of initial conditions is applicable, and succeeding 
positions indicate the particular parameters that are 
desired in each case. Ordering of parameter entries 
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for IC Types 1, 2, and 3 is as shox«m in Section 
3.3,2. 

b. DPRAM (Differential Equation Parameters) 
(Lines 16 through 20) 

The ordering of DPRAM differential equation para- 
meter entries is as shox«m in Section 3,3.2 • How- 
ever, in this application the Ti and T2 para- 
meters are associated with Satellite 1 only, 

c. PSAT2 through PSAT6 (Initial Condition and Drag 
Coefficient Parameters for Satellites 2-6) 
(Lines 21 through 25) 

The ordering of PSAT IC and drag parameter entries 
for Satellites 2 through 6 is shoxm below. 



PSAT2 


^2 


62 


^2 


A2 


^2 


V2 


t02 


DRAG2 


PSAT 3 


^3 


^3 


33 


A3 


r3 


V3 


t03 


DRAG 3 


PS AT A 


Ch 


64 


ei* 


A^ 


rit 


vi+ 


to. 


DRAG, 


PSAT5 


^5 


«5 


es 


A5 


rs 


V5 


tos 


DRAG 5 


PSAT6 


^S 


«6 


ee 


As 


^6 


V6 


^H 


DRAG 6 



It is important to note that only IC Type 2 may be specified for 

Satellites 2 through 6, 

d, KPRAM (Orbit Parameters for Satellite 1) 
(Lines 26 and 27) 

Ordering of KPRAM orbit adjust parameters is as 
previously shox«m in Section 3,3,2, It should be 
noted that the orbit adjusts are for Satellite 1 
only. 

Sixty trajectory parameters is the maximum number 
which may be selected for simultaneous solution. 
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3.5.2.2 Radar Parameters 



28 


M 


RAPAR 


05, 99 


































29 


D 


01, 01 


B 


N 






L 


A 


T 












30 




03, 01 


.01 




31 




04, 01 
















32 


D 


01, 02 


C 


K 


2 


4 


R 


B 


I 


A 


S 








33 




03, 02 


500. 




34 




04, 02 


750. 














35 


D 


01, 03 


H 


U 






T 


B 


I 


A 


S 


- 






36 




03, 03 


.001 




37 




04, 03 


.0035 





Line 28 indicates that the data listed subsequent to RAPAR on the 
load sheet will be input into a 5 x 99 matrix array which has been 
preset to zero. The columns of this array correspond to parameters, 
and the rows correspond to parameter identification (Positions 1 and 
2), bounds (Position 3), and bias estimates (Position 4) respectively. 
Row 5 is currently not used. 

Lines 29 through 31 specify the first applicable radar 
parameter. Line 29 contains the station name (Positions 1 and 2), 
pass identification (Positions 3 and 4), and parameter name (Positions 
5 through 10). Lines 30 and 31 contain the bound and initial value, 
respectively, in feet, degrees, and minutes except for R, P, and Q, 
which are in feet per second. If the parameter is station latitude, 
longitude, or altitude, the initial value is taken from the station 
location card. 
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If the pass identification character position is left 
blank, all data with the indicated station name will be used to cor- 
rect the parameter. If the pass identification is not omitted, only 
data that are identified by both the indicated station name and indi- 
cated pass identification will be used to correct the parameter. If 
the radar parameter specified is station latitude, longitude, or alti- 
tude, the pass identification is ignored and all data with the station 
name are used for the parameter correction. 

Lines 32 through 34 and Lines 35 through 37 specify the 
second and third radar parameters respectively. Station names, pass 
identifications, parameter names, bounds, and initial values are 
treate^ in the same manner as the inputs in Lines 29 through 31 des- 
cribed above. Additional cards may be added in cases where more than 
three radar parameters are involved. Available radar parameters are 
listed in Table 9, 

Note that the total number of parameters which may be se- 
lected for simultaneous solution must be less than one hundred. 
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Table 9. 



PARAMETER 


SYMBOL 


Station latitude 


LAT 


Station longitude 


LONG 


Station altitude 


ALT 


Time bias 


TBIAS 


Range bias 


RBIAS 


Azimuth bias 


ABIAS 


Elevation bias 


EBIAS 


Topocentric right ascension bias 


RTBIAS 


Topocentric declination bias 


DTBIAS 


Topocentric hour angle bias 


HABIAS 


Geocentric right ascension bias 


RGBIAS 


Geocentric declination bias 


DGBIAS 


Argument of latitude (u) bias 


UBIAS 


Cross plane (v) bias 


VBIAS 


Height bias 


HBIAS 


X bias 


XBIAS 


y bias 


YBIAS 


z bias 


ZBIAS 


P bias 


PBIAS 


Q bias 


QBIAS 



Range-rate bias 

P bias 

Q bias 

Range scale factor (Kr) 

Range-rate scale factor (Kd) 



RDBIAS 

PDBIAS 

QDBIAS 

KR 

KD 
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3.5.3 Lines 38-50 Parameter Bounds 



38 




ENDS 


.5 




39 




2 


.5 




40 




3 


.1 




41 




4 


.5 




42 




5 


1000. 




43 




6 


5 




44- 




7 


.05 




45 




8 






46 




9 






47 




10 






48 




11 






49 




12 






50 




13 







A bound must be entered for each parameter selected. These 
bound entries must be in the same sequence as the parameters. For 
each iteration of the differential correction process, the change in 
each parameter is less in absolute value than the corresponding bound 
if that bound is positive, zero if the corresponding bound is zero, 
or unrestricted if the corresponding bound is negative. 
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3. 5. A Lines 51-62 Observation Slgmas 



51 




SIGMA 


100. 




52 




2 


.5 




53 




3 


.5 




5A 




A 


1000. 




55 




5 


200. 




56 




6 


300. 




57 


I 


ISIG 


1 




58 


I 


2 


2 




59 


I 


3 


3 




60 


I 


A 


113 




61 


I 


5 


llA 




62 


I 


6 


115 





The observation sigmas for each observation type must be input 
in the same manner as was previously well described in Section 3.4.3. 

3.5.5 Line 63 Maximum Iterations 



63 


I 


MAXIT 


A 





The symbol MAXIT must be input with an Integer value' to specify 
a stopping point to the run. If the differential correction process 
has not converged at the end of MAXIT iterations, the process will be 
terminated. 

3.5.6 Lines 6A-66 Residual Editing 



6A 


CEDIT 


3 




65 


2 


.9 




66 


3 


1 





61 



--T^ 



If CEBIT is zero, no editing is done. If CEBIT is non-zero, 
residuals will be edited in accordance with the following: 

a. CEBIT ^ Residuals greater than (input 

Sigma X I CEBIT I ) will be discarded. 

b. CEBIT -^ Residuals greater than (statistical 

sigma from the previous iteration >^ CEBIT) 
will be discarded. No editing is done on 
Iteration 1. Sigmas are computed for the 
first five data types encountered for each 
station. 

Line 65 represents a scale factor such that if CEBIT (2) is non- 
zero, CEBIT is replaced by (CEBIT x CEBIT(2)) at the end of each itera- 
tion. If CEBIT(2) is zero, CEBIT is not modified. 

Line 66 is a special option wherein if CEBIT (3) is non-zero. 
Iteration 1 will be repeated with editing performed with the sigmas 
computed during the first pass through Iteration 1. This will allow 
editing to be done on all iterations with computed sigmas. 

Sigmas (rms) for the first five data types for each station are 
computed and printed at the end of each iteration regardless of the 
residual-editing option selected. 

3.5.7 Bata Correction Factors 

3.5.7.1 Lines 67 and 68 Elevation Angle Refraction 



67 


1 


REFR 






, 

! 68 


1 1 ■! 


2 


4.0 


-6 



A table of refraction indices Hi, which may contain 
up to ten values, may be input starting at REFR. The entry used to 
compute refraction corrections for radar elevation observations is de- 
termined by the type number contained in Column 6 of the corresponding 
station location card. A zero in Column 6 causes the entry at REFR to 
be used, a 1 in Column 6 causes the entry at REFER -f- 1 to be used, etc. 

If the table contains no entries, the value 
312.0 ^ 10"^, which is built in at location REFR, will be used to 
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compute refraction corrections for all data whose station location 
cards contain zero in Column 6. All other positions of the table are 
assembled as zeros. 

3.5.7.2 Range Refraction Correction 



69 




RREFC 


1 





Refraction corrections to all range observations will 
be computed and applied if RREFC is non-zero. 

3.5.7.3 Line 70 Propagation Time Correction 



70 




SLT 


-2820.1763 





The velocity to be used in calculating the observation 
time correction due to propagation time is entered at SLT in units of 
earth radii per minute. In the absence of an entry, no correction will 
be applied. If an entry is present the correction will be applied to 
times associated with R, A, E, R, h, P, Q, P and Q data only. 

3.5.8 Lines 71, 72 Data Tape Options 



71 


I 


IBCDI 


6 




72 


I 


IBINI 


1 





If the radar observation and station location information is to 
be input via a BCD tape other than the normal FORTRAN system input tape, 
the tape number must be specified at IBCDI as 11. If a binary tape 
containing compacted radar data produced by a previous run is to be 
input, IBTNI must be specified. 
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3.5.9 Lines 73, 74 True Equinox Correction 



73 




DALFG 


-.004 




74 




2 

















An additive factor may be applied to the computed right ascen- 
sion of Greenwich at midnight of epoch day by entering the appropriate 
value in units of degrees at DALFG for Vehicle 1, at DAFLG(2) for 
Vehicle 2, etc. This entry usually is used to correct from mean to 
true equinox reference coordinates. 

3.5.10 Lines 75-78 Proximity Indicator 



t* 75 




AN0M1 




A 76 




2 




h 77 




3 




Rq 78 




4 





During an orbit determination run, an Indicator may be obtained 
whenever the trajectory passes within a given distance (range) of a 
point on the surface of the earth by input of geodetic latitude (deg) , 
east longitude (deg), and altitude (n mi) of the point at ANOMl and 
the succeeding two positions and of the testing distance (the range 
from the point to the vehicle) at AM0N1(4). Testing and printing is 
done for up to three such sets of input at AMONl, AN0M2, and AN0M3. 

3.5.11 Lines 79 on Parameter Constraint Options 

There are two different parameter constraint options available 
in TRACE-D. Both are described below in detail. 

3.5.11.1 General Constraint Matrix 

Option 1 allows the user to apply linear constraints 
to any of the variables in the parameter set, and to specify these 
constraints by means of an input matrix set up as shown below. 
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79 


I 


KNST 4 




80 


I 


BLIST 1 




81 




1 




bli 82 




1 




83 




2 




84 


I 


2 




b?? 85 




1 




86 


I 


3 


87 




3 ! 


^3? 88 




1 


89 




4 


90 




4 


b,, 91 




1 

i 


92 


, 


5 1 

> 1 


93 




! 


bsi 94 






95 




6 ( 


96 




2 


be?. 97 




•5 ! 


98 




5 




99 


I 


5 




C51 100 


1 


6 




101 


1 ^■ 


7 




102 


I 


5 




103 




i 



The clearest explanation of how this constraint matrix 
works is an example. Therefore, in the case shown above it is as- 
sumed that n parameters are to be solved for where 
(pi, P2, ..•» p^) = P* T^^e ordering of the Pi corresponds to the 
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order of the X's for the CPRAM, DPRAM, and KPRAM and the RAPAR ar- 
rays, respectively. Further assuming that these parameters are to be 
subjected to m linear constraints, which, for example for n « 6 
and m = 2 might be Pi + P5 = 6, P2 - 2p5 « 0, KNST above would be 
equal to (n - m) = 4, or the number of effective (unconstrained) 
parameters. 

The BLIST constraint matrix is represented by the factor 
B in the expression p « Bp + c , where the p are the effective 
parameters. In connection with the foregoing example, this expression 
would assume the form 



Pl 




P2 




P3 


a 


Pi+ 




P5 




P6 





1 














1 














1 














1 


- 1 














.5 








~v 






->' 





r ~ 
Pl 









P2 







X 


P3 


+ 







Pk 





6 





L -1 




^^ ^ 



V. 



B 



c 



Pertinent constraints are applied by input of the augmented (n + 1) 
by (n - m + 1) matrix 

B 





:] 



in the BLIST array. If aij is an element of this augmented matrix, 
then i, j, and a^-j would be input. The inputs for the example 
given above would be as indicated in Lines 79-103. KNST (line 70) is 
input as 4 and then the BLIST matrix follows. Lines 80, 81, and 82 
specify that BLIST (1, 1) is a 1 and lines 83, 84, and 85 that BLIST 
(2, 2) is a 1, etc. Elements of the matrix not input are assumed to 
be zero. In this way any subset of the set of parameters may be 
linearly constrained and the constraint matrix may be inout to the 
program in the exact same manner. 
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3.5.11,2 Special Application to Station Location 

The other option available is a specific case of the 
general constraint matrix which is applicable in the following type of 
problem. Consider a differential correction run with data from a tri- 
static system and a set of parameters including the latitudes, longi- 
tudes and altitudes of the three stations. In such a case, by choosing 
this constraint option^the correction of the station coordinates can 
be performed such that the baseline distances are held constant. 

Mathematically, this implies that three constraining 
equations^ are applied to the nine station locations; more specifi- 
cally they are of the following form. 

A4)i = ai Aei + a2 AHi + as Ae2 + ai4 AH2 + as A(f)3 + ag AH3 

A(j)2 = bi Aei + b2 AHi + b3 Ae2 + bi+ AH2 + bs A4)3 + bg AH3 

Ae3 = ci ASi + C2 AHi + C3 Ae2 + Ci+ AH2 + C5 A(f)3 + cg AH3 



where 



(ei, *i. Hi) 



are the three sets of latitude. 



,Q . „ >. Longitude, and altitude for 
(92, *2, H2) Stations 1, 2, and 3. 

(63, 4*3, H3) 

and a., b^, c. i = 1, ..., 6 are the 18 coefficients of the 

three equations. 

With the equations in this particular form, it seems 
that the general constraint option could be applied easily. However, 
there are two reasons for the separate, treatment of this application 
that are not immediately apparent. First, the constants a^, b-j^, c^ 
i = 1, ..., 6 are functions of the current station location^ and 
are not easily and accurately calculable by the user. Second, since 



1 

The detailed derivation of this particular form of the constraining 

equations is not given here but may be found in Reference 3. 
2 

The exact expressions for the coefficients can be found in Reference 3, 
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the constants are functions of the current station locations, they 
must be reevaluated every time the station locations change (i.e., 
after every iteration of the correction process the constants must be 
reset since the station locations have been adjusted). In other words, 
the type of constraint matrix needed here is a dynamic one that is 
adaptable to change throughout the estimation procedure. It is now 
obvious that the general constraint matrix of Section 3.5.11.1 is too 
rigid for this purpose. 

Thus, added to the TRACE-D program is the capability to 
compute these necessary constants from input information, insert them 
into the appropriate location of the constraint matrix, recompute the 
constr *"s as a function of the new station locations after each itera- 
tion an J adjust the constraint matrix accordingly. Once the user has 
set up the appropriate input the rest is performed by the program. 

The input necessary is best described with an example. 
Consider a differential correction run with data from three stations 
and a parameter set consisting of the ADBARV elements, drag, and sta- 
tion locations. To employ the fixed baseline distance constraints on 
the station locations consider B in the following expression. 

Ap' = B AP + c 

There are 16 parameters to be estimated but there are 3 constraint 
equations, therefore there are 13 effective or unconstrainted para- 
meters. B is then of dimension 16 x 13 and looks like this: 
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Ao 
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M 




AE 




A/ 


L 


AE 




AV 
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A 


drag 


A 


latj 


A 


longi 


A 


alti 


A 


lat2 


A 


long2 


A 


alt2 


A 


lata 


A 
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A 


alt3 
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ai 


32 


33 


34 


as 


ae 

















1 



































1 





























bl 


b2 


b3 


b. 


bs 


be 























1 



































1 



































1 




















ci 


C2 


C3 


C4 


C5 


C6 




















_0_ 











_1_ 



[Ap] 




[c] 


Aa 







A6 







AB 







AA 







AR 







AV 







A drag 







A longi 


+ 





A alti 







A long2 







A alt2 







A lat3 







_A_alt3 

























The quantities a^, b^, Cj^, i = 1, ..., 6 are, of course, the coef- 
ficients in the constraint equations which are unknown to the user. 

This B matrix is input in exactly the same manner as the 
general constraint matrix except that each element that is a coeffi- 
cient (either ai, b^^, c-^) must be input with a value of - 10,000. 
This is merely to indicate to the program the position in the matrix 
of these terms. 

Lines 104 on show the appropriate format for the necessary 
input of this case. 
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104 




NUMB 




105 ] 




44 


-1 


106 ] 




KNST 


13 


107 ] 




BLIST 


1 


108 ] 






1 


109 






1 


110 ] 


[ 




2 


111 ] 


[ 




2 


112 : 






1 


113 








114 ( 


itC. 






115 









116 




8 




117 




8 




118 




-10,000. 




119 




8 




120 




9 




121 


etc. 


-10,000. 





Lines 104 and 105 specify the quantity NUMB (44) which must be input 
as a - 1 to choose this constraint option. Line 106 is the number of 
unconstrained parameters in the parameter set. Lines 107 on are the 
non-zero elements of the B matrix entered in order. 

3.5.12 Station Identification Cards 

The locations of tracking stations or of points on the surface 
of the earth which are associated with observations must be input by 
means of station cards carrying appropriate information indicating 
the manner in which the corresponding observations are to be processed. 
These cards were previously discussed in Section 3.4.5 and shown in 
Figure 1, but the format is included again here for completeness. 

a. Columns 1 and 2 (ST) : Station identification symbol. 
No two stations may be identified by the same symbol 
or any one station by the symbol TS. 

b. Column 5: Sigma index identifying observation-sigma 
set to be applied to data from corresponding station. 
The sets of sigmas are input with the base data. 

c. Column 6: Type of refractivity correction to be used 
for elevation readings from this station. Refracti- 
vities are numbered in their input order within the 
base data (see Line 67, REFR) . 
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d. Columns 9 througli 17 : North latitude of station in 
degrees. 

e. Columns 1 9 through 27 : East longitude of station in 
degrees. 

f. Columns 29 through 36 : Altitude of station in feet. 

g. Columns 38/39 and 41/42 : If a station reports P, 0, 
or f*, data, Columns 38/39 and 41/42 contain the 
two letter symbols for the associated station (s) of 
the tracking configuration. Each such associated 
station must be represented by a separate station 
card, but it is not necessary for Columns 38/39 
and 41/42 to be filled out on the latter. 

These cards must be preceded by a card with END in Columns 1-3 and 

followed by a card with TS in Columns 1 and 2. 

3.5.13 Observation Data Cards 

The observations to be processed must be the last set of cards 
in the deck and must follow the station location cards. 

Since the number of observations that may be input to TRACE-D 
is unlimited, there is a flocking procedure that must be followed. A 
flock of data consists of 200 or less observation cards for a single 
satellite. There is no limit to the number of flocks that can be 
stacked; however, each flock must be followed by a card with TF in 
Columns 1 and 2 and the times associated with the observations of a 
flock can never be earlier than the latest time of the previous flock. 
This means that the data is not time sorted within a flock but the 
flocks are sorted according to earliest and latest times therein 
contained . 

The contents of an observation card is formatted in the fol- 
lowing way, as shown in Figure 4. 

a. Columns 1 and 2 (ST): Station Identification, symbol 
which must correspond to a station location card. 

b. Columns 3 and 4: Pass identification (optional). 

c. Columns 5 through 21: month (columns 5 and 6), day 
(columns 7 and 8), hours (columns 9 and 10) 
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minutes (columns 11 and 12), and seconds (columns 
13-21) , of the corresponding observations relative 
to Greenwich Mean Time. 

d. Column 22: Observation set number. 

e. Columns 23-76: Observation 1 (columns 23-40), 
Observation 2 (columns 41-58) , Observation 3 
(columns 59-76) . 

f. Column 77: .Card number indicating observations, 
variances, or covariances. 

The observation set (d) number and the card number (f) are the two in- 
dicators that define exactly what type of observations are contained 
on each card. Table 10 gives a complete listing of all types according 
to obs'^rvation set number and card number. 

After the last observation card in the deck, a TF card is not 
necessary, however a card with TR in Columns 1 and 2 plus an END card 
must be attached as the last two cards of the deck. If data from more 
than one satellite are to be used, the Satellite 1 data are set up in 
the same manner except that the TR card must be replaced by a TT card. 
The data for Satellite 2 are then similarly arranged. If Satellite 2 
is the last vehicle, corresponding data are followed by a TR card and 
an END card: if it is not the last satellite, data are followed by a 
TT card. Data for Satellites 3 through 6, as applicable, are added in 
the same manner. A TR card rather than a TT card must follow the data 
corresponding to the last satellite, and an END card must follow the 
TR card. See Section 3.6 for a picture of the correct setup. 

3. 6 Input Deck Arrangements 

Since the basic running deck to be used for each function to 
be performed can be obtained, the user need not concern him- 
self with its contents. Therefore in each of the Figures the separate 
parts of this deck are ignored and it is labeled totally as the basic 
running deck. The categories of data used are described in detail in 
previous sections (for instance, Basic Data is explained in Section 
3.2, etc.). 

72 



CO 

.-I 

^4 



0) 
•H 






c 




/"-N 


o 




r^ 


•H 


(U 


r^ 


4-1 


p 




m 


>. 


C2 


> 


H 


g 


U 




D 


0) 


T^ 


.H 


03 


H 


O 


JO 


m 


U 


o 


U 


'^.^ 


c 




x^s 


o 


u 


CN 


•H 


Q) 


CN 


*J 


JO 




m 


e 


c 


> 


D 


g 


H 


2 


D 


(U 




.H 


03 


4J 


O 


XI 


(U 


U 


O 


C/3 


v^ 









^-^ 
































/"-N 


? 


^ 














^„^ 




/-^N 






/-^N 




/'~\ /'~\ 




w 


#< 














X 




< N 




^-^ 


O" 




cy -o- 




#< 


'>*-• 


H 














»> 




#< 




cr 






"*-• #> 


/"-N 


<c 


/-^N 


<o 














> 




<>-. 




K^ 


Qm 




W -PL, 


W 


•-N ^— ' 


O Q) ^ 


""•^^ 










^•-\ 


/"-N 


"*-• 


^•-\ 


'^.^ 




Q) 


^ ^— ' 




(U •-N ^-^ 


'>*-• 


W (U 


^ r-{ :ii 


Q) 










x: 


x: 


OJ 


< N 


(U 




O 


o- (U 


Q) 


O'Cy (U 


c 


^-^ o 


J-I tJt'^-' 


O 










'^.^ 


'^.^ 


o 


'^.^ 


o 




c 


^-^ o 


4J 


c ^-^ o 


o 


(U C 


4-1 C (U 


C 










OJ 


OJ 


c 


(U 


c 




Q) 


(U c 


m 


(U (U c 


•H 


o m 


c m o 


m 










-o 


o 


m 


o 


m 




U 


o m 


J-I 


J-I O CO 


U 


C -H 


(U c 


*H 










D 


c 


•H 


c 


•H 




Q) 


C ^H 




(U C -H 


m 


m J-i 


o v-i m 


V-i 










4-1 


m 


J-I 


m 


J-I 


OJ 


M-I 


m J-I 


(U 


M-I rt J-I 


> 


*H m 


O D *H 


m 










'H 


*H 


m 


-H 


m 


6;. M-i 


•H m 


60 M-< -H CC 


(U 


J-i > 


a o n 


> 










4-) 


J-I 


> 


J-I 


> 


C 


*H 


J-< > 


C 


•H J-I > 


iH 


m O 


o j:: m 


O 










iH 


m 


O 


m 


O 


m 


-o 


m O 


m 


-O CO O 


P^ 


> U 


H > 


u 










1 


> 


U 


< N > 


U 


Ci^ 




> u 


Ci^ 


> u 




"""I 


^^ 






/"-N 






r 






















/-^N 


H 






6£ 










/— s 




/-^N 






^-\ 




/-^N 




w 


<o 






to 






/"-N 




J-I 




C N 




y-^ 


O' 




^ -o- 




#< 


v*^ 






s^* 






> 




#< 




#< 




P^ 


*^ 




.(X, 




Ci^ 


C ^N 






C 


/-^N 




K^ 




D 




X 




'^.^ 


^ 




w .fii 


/"-N 


•-N ^^ 


o o E- 






o 


6C 


Q) 


/"-N 


'>*-• 


/"-N 


"*-• 




(U 


^ ^— ' 




(U ^-N^^ 


<c 


< (U 


-1-i '1-{ <<i 




O 


•H 


K3 




C 


> 


(U 


< >. 


(U 




o 


(X, (U 


(U 


O -(X, (U 


"*-• 


v-' o 


u u ^^ 




•H 


4J 


"*-• 




m 


"*-• 


o 


^-^' 


o 




c 


v-' o 


4-J 


C ^^ o 




(U c 


4J m (U 




V-i 


m 


0) 




iH 


(U 


c 


(U 


c 




(U 


(U c 


m 


(U (U c 


j:: 


o m 


ceo 




4J 


c 


o 




a 


o 


m 


o 


m 




J-I 


o m 


J-I 


u o m 


4-J 


C -H 


(U -H C 




c 


•H 


c 






c 


•H 


c 


•H 




(U 


C -H 




(U C -H 


D 


m U 


o ^ ^ 




(U 


.-I 


m 




03 


m 


J-I 


m 


J-I 


(U 


M-< 


m J-I 


(U 


M-I m J-I 


e 


•H m 


O O -H 




o 


O 


*H 




tn 


•H 


m 


•H 


m 


60 M-I 


H m 


60 M-I -H m 


•H 


H > 


O (U J'Ht 




o 


(U 


V-i 




o 


J-I 


> 


J-I 


> 


C 


•H 


J-I > 


C 


•H J-I > 


N 


m O 


o xi m 




(U 


-o 


m 




J-I 


m 


o 


m 


O 


CO 


xt 


m O 


m 


-o m O 


<c 


> u 


H > 




o 




> 




u 


> 


u 


< >^> 


U 


Ci^ 




> u 


Ci^ 


> u 






AJ 


^,_^ 








/"-N 


























j:: 


H 


4-1 






6C 
























/-^N 


60 


K3 


j:: 






KD 






/"-N 




/"-N 






/"—^ 




/"-N 


/-^N 


<c 


•H ^N 


»> 


60 


/"-N 




#« 






> 


< 


>% 


^-s 




(X, 




•(X, 


Pti 


#< 


H H 


H 


•H 


6C 




60 






#< 




»v 


Pi 






/— s 




'>*-• 


Ci^ 


a ^ 


C5 


u 


a 


/"-N 


a 


/-^N 




D 


< 


X 


"*-• 




a: 


•Ci^ 


•Ci^ 


(U 


•-N ^-^ 


o ^^ H 


'>*-• 






6C 


^-^ 


M-i D 


/-^N 


'>*-• 


/— s 


'>*-• 


(U 




/"-N "'-^ 


'>*-• 


/-^N ^-^' 


60 Di; (U 


•H C C5 


(U 


o 


c 


a 


OJ 


O ^-^ 


D 


(U 


< X 


(U 


60 




^ OJ 


(U 


•P^ (U 


c 


^^ o 


V-i o ^-^ 


o 


•H 


o 


'>*-• 


o 


(U 


'^.^ 


o 


'>*-• 


o 


C 




v-' o 


4J 


v-' o 


m 


Q) C 


4J -H (U 


c 


V-i 


♦H 


OJ 


c 


W -o 


(U 


c 


(U 


c 


en 




(U c 


m 


(U c 


n 


O OJ 


c w o 


c« 


4J 


W 


o 


m 


c a 


o 


m 


o 


m 


J-I 




o m 


J-I 


o m 




C -H 


OJ c c 


•H 


c 


C 


C 


•H 


0) 4-) 


c: 


»H 


c 


•H 






C -H 




C -H 


4-1 


m V-i 


O (U m 


H 


(U 


Q) 


m 


V-i 


E ^ 


m 


J-I 


m 


J-I 


4J 




m J-I 


0) 


m J-I 


c 


•H m 


O O 'H 


m 


o 


O 


*H 


m 


a w 


*H 


m 


•H 


m 


c 




•H m 


60 


•H m 


m 


M > 


CL W J-l 


> 


o 


W 


V-i 


> 


60 m 


J-I 


> 


J-I 


> 


m 




J-I > 


c 


J-I > 


.H 


m O 


O CO m 


O 


(U 


c« 


m 


O 


J-i t-H 


m 


O 


m 


O 


.-I 




m O 


m 


m O 


en 


> u 


H > 


u 


a 




> 


u 


< 


> 


u 


< X > 


u 


CO 




> U 


p^ 


> u 


o 


.H CM 


O i-H 


CN 


o 




.-I 


CN 


O 


.-I 


CM 


O i-H 


Cvl 
1 
1 


o 




.-I CM 


o 


.-I Cvl 


.-I 


.H .H 


CN CN 


CN 


CO 




CO 


CO 


<r 


<r 


<r 


in m 


in 


vO 




\C \D 


r^ 


r^ r^ 



03 

J-I 
O 

o 

c 



c 

03 

(U 
D 

H 

C 

o 

•H 

> 
U 
(U 
03 
X 
O 



X 
H 



73 



3.6,1 Trajectory Prediction Deck Arrangement 

When the function to be performed is a single trajectory pre- 
diction (function code ITIN = 3), the arrangement is as in Figure 5. 



c 



END 



END 



TRAJECTORY ONLY DATA 



BASIC DATA (CONSTANTS, 
INITIAL CONDITIONS, DYNAMIC 
MODEL etc.) 



T SUBTYPE, DATA 



BASIC RUNNING DECK 




Figure 5. 
TRAJECTORY PREDICTION DECK ARRANGEMENT 

When several trajectory predictions are to be done in one run 
(function code ITIN = 333, etc.) then the deck arrangement is merely 
an extension of that for .the single case (Figure 6) • 

3.6,2 Observation Generation Deck Arrangement 

Figure 7 and Figure 8 show the correct data setup for a single 
data generation run (function code ITIN = 4) and several cases (func- 
tion code ITIN ^ 4.4, etc.) respectively. 
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Figure 6. 
MULTIPLE TRAJECTORY PREDICTIONS DECK ARRANGEMENT 
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Figure 7. 

DECK ARRANGEMENT FOR 
OBSERVATION GENERATION 
SINGLE CASE. 



TR 



DATA TYPE CARDS 



TR 



\ 




DATA LIMIT CARDS 



/ 



^ 



TS 



STATION CARDS 



^ 



END 



END 



OBSERVATION GENERATION DATA 



BASIC DATA 



T SUBTYPE, DATA 



BASIC RUNNING DECK 




J/ 



I 



J 
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TR 



Figure 8. 

Deck Arrangement for 
Two Cases of Observa- 
tions Generation 



DATA TYPE CARDS 



[ 



TR 



DATA LIMIT CARDS 



[ 



TS 



STATION CARDS 



END 



^ 



^ 



y 



SECOND SET OBS. GEN. DAT/ 



[ 



TR 



DATA TYPE CARDS 



J. I 
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L 



TR 



DATA LIMITS CARDS 



^ 



TS 



r 



STATION CARDS 



END 



END 



OBSERVATION GEN. DATA 



/ 



/ 



h 



BASIC DATA 



T SUBTYPE, DATA 



\ 



BASIC RUNNING DECK 



y 



y 



J 



]/ 



J/ 
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3.6.3 Orbit Determination Deck Arrangement 

There are several different deck setups that are pemrlssible 
for an orbit determination run (function code ITIN = 12) according to 
the choice of several suboptions. 

Figure 9 is the simplest deck setup and applies when the ob- 
servation data is being input on tape (i.e., IBCDI input variable is 
used) . 

Figure 10 is the ordinary deck setup for a single case when 
the observation data is attached to the deck. Figure 11 is the deck 
setup for two such stacked cases. 

Figure 12 is the appropriate data arrangement for the situa- 
tion of multiple satellites (in this case three satellites). 




Figure 9. 

Deck Arrangement for Orbit Determination with Input Tape* 

Note that the BCD tape of observations must contain as its first 
set of cards the station cards followed by a TS card and then the 
observations in proper flocks followed by TF cards with a TR card 
replacing the TF card cf the last flock. 
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Figure 10. 

Deck Arrangement for Single 
Orbit Determination Run with 
All Data on Cards 
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END 
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TF 
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END 
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*/ TRACKING ONLY DATA 
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T SUBTYPE. DATA 
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/ / 
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./ 
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* Note that the word "tracking' is used interchangeably with the words 
"Orbit Determination". 
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Figure 11. 

Deck Arrangement for 
Two Orbit Determina 
Cases All Data on 
Cards 
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Figure 12. 

Deck Arrangement for Multiple 
Satellite Orbit Determination 
Run (3 Satellites) 
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3-6.4 Sample Input Decks 

Figures 13, 14, and 15 are listings of sample data decks for 
trajectory prediction, observation generation, and orbit determination. 
Besides giving the user an idea of what a typical deck looks like, 
they may also be used to familiarize the user with the operation of 
TRACE-D before he creates his own input. Another reason for running 
a test case of this sort is to determine the status of the program 
tape, for it is possible for the contents of a tape to deterioriate 
and no longer be useful . 

It is not to be assumed that these sample cases exercise all 
the possible options of TRACE-D, but rather that they exhibit typical 
input -^d output. A complete description of the output obtained from 
these runs is included in Section 4. 



^ In this instance, there is no cause for alarm. MITRE has the capa- 
bility to create a new and correct program tape in 24 hours. 
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Section 4. 

SAMPLE OUTPUT 
4.0 This section presents the actual output from the test cases 
described in Section 3.6.4 whose input is listed in Figures 13, 14, 
and 15. Each test case is covered separately and the code numbers 
next to certain sections of the ouotut correspond to the numbered 
paragraphs of the explanation that follows . 

Due to lack of space only pieces of the total output are 
shown; however, the full actual output is available. 

^ • ^ Trajectory Prediction Output 

The following sequenced paragraphs describe the different 
parts of trajectory prediction test case output shown in Figures 16 
to 19. 

1. The constants used in the program (c, INTEG, NUMB, 
IFLAG) are printed out here in exact input card 
image . 

2. These numbers are the geopotential inputs printed 
out exactly as read in. 

3. The last part of the input data printed out is all 
the traj ectory-only input . 

4. Program name and identification of the particular 
Aerospace version. 

5. The Epoch time for this run. 

6. Indicators printed out givine the logical flow 
from segment to segment and old and new func- 
tion numbers. 

7. Card image of last END card in the deck. 

8. Trajectory initial conditions in three coordi- 
nate systems. Initial-condition values as 
shown are the result of transformations which 
have been applied to the input values. The 
transformation for the input coordinate set 
(a. 6, 3, A. r, v in this case) consists of 
conversion from decimal to octal numbers, 
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conversion of units from feet, degrees, and 
seconds to earth radii, radians, and minutes, 
and performance of corresoondine inverse non- 
versions for output. The two other tynes of 
elements sets also require accomplishment of 
coordinate-system transformations in addition 
to the number- and units-systems conversions 
noted above. Accuracy of the values as 
printed therefore is subject to numerical 
roundoff errors. 

Ouantities in the left-hand column are 
position and velocity components in the basic 
vernal equinox coordinate system, with units 
of feet and of feet oer second. The center 
column gives the usual ADBARV spherical sys- 
tem coordinates (i.e., Type-2 initial condi- 
tions in units of feet, feet per second, and 
degrees) . The right-hand column from top to 
bottom contains orbit semi-major axis, eccen- 
tricity, inclination, right ascension of as- 
cending node, argument of perigee, and time 
of last perigee passage in minutes from mid- 
night of epoch day. Other units are feet and 
degrees. 

9- Identification of atmosphere model to be used 
in computing drag force. LOCKHEED and ARDC 
1959 refer to the Lockheed-Jacchia and the 
ARDC 1959 model atmospheres, respectively. 
Reciprocal of ballistic coefficient, CdA/W. 
Ballistic coefficient , W/Cj^A. dj and d2 
are values of certain constants in the 
Lockheed-Jacchia atmosohere density expres- 
sions. 

10. The product GM, or the universal gravita- 
tional constant times the mass of the earth 
(frequently designatf^.d by y) , expressed in 
units of earth radii cubed per minute sauared- 
Unitless coefficients of the zonal harmonic 
terms in the earth potential field expansion 
(i.e., J 2 through Ji 5) • 

Coefficients and longitudinal arguments of 
tesseral and sectorial terms in earth poten- 
tial field expansion, with J indicating a 
coefficient and L an argument in desirees. 
The digits following J or L are the 
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associated Legendre polynomia!!. degree and order, 
respectively. 

11. Names of solar-system bodies included in computa- 
tion of perturbative accelerations. In this case 
none are used. 

12. Trajectory method indicator. The Cowell formula- 
tion of the equations of motion is the only tra- 
jectory method available in the present TRACE-D 
program. 

Numerical integration parameters. The Gauss- 
Jackson method (subroutine DE02) is the only 
integration method available in the present 
TRACE-D program. 

13. The indices Ni and N2 input in the lOB 
array respectively indicate the highest-degree 
zonal term and the highest-degree and -order 
tesseral term. No term higher than those indi- 
cated by these input numbers will be computed, 
regardless of the model coefficients and argu- 
ments that are input to the physical constants 
region. 

14. Physical constants and units conversion factors. 
Decimal and octal equivalents are given for the 
following quantities. 

a. OMEGA E(oJe) Earth rotation rate in 

radians per minute. 



b. ALPHA G(ag) 

c. FT/KM 

d. E.R/A.U. 

e. I-O, DISTANCE 

f . G 



DEG/SEC// 
RAD/MIN 



Right ascension of Greenwich 
at midnight of the day of 
epoch in radians. 

Feet per kilometer. 

Earth radii per astronom- 
ical unit. 

Input/output units conver- 
sion factor in feet per 
earth radius. 

Acceleration due to gravity 
in feet per second squared. 

Units conversion factor in 
degrees per second/radians 
per minute. 
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h. RELATIVE 

MASS (SUN) 

i (VENUS) 
j . (JUPITER) 
k. GM 

1 . EARTH 

RADIUS-FT 

m. N.M/E.R 

n. I-O, VELOCITY 

o. FT/SEC// 
E.R./MIN. 

p. 1/EPS 

q. (MOON) 

r. (MARS) 

s. (SATURN) 



Mass of sun relative to 
mass of earth, 

Mass of Venus relative to 
mass of earth, 

Mass of Jupiter relative 
to mass of earth. 

Earth gravitational constant 
(u) in earth radii cubed per 
minute squared. 

Number of feet per earth radius 
for internal calculations 
(rather than I-O) , 

Number of natuical miles per 
earth radius. 

Conversion factor for input/ 
output in feet per second/ 
earth radii per minute. 

Units conversion factor in 
feet per second/earth radii 
per minute. 

Earth flattening (reciprocal 
of earth ellipticity) , 

Mass of moon relative to 
mass of earth, 

Mass of Mars relative to 
mass of earth. 

Mass of Saturn relative to 
mass of earth. 



15. Date and time of day (Greenwich Mean Time) with which the 
quantities in the print block following are to be asso- 
ciated. 

16. Minutes from epoch, minutes from midnight of current day, 
system time (i.e.^ seconds from midnight of current day), 
and current integration step size in minutes. 

17. X, y, z, r: Componenents and magnitude of the radius 
vector from geocenter to satellite in the basic coordi- 
nate system in units of feet. 

18. X, y, z, v: Components and magnitude of the inertial 
velocity vector with respect to the basic coordinate 
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system in units of feet per second. 

19. Geodetic latitude, in degrees, of the point where the 
radius vector intersects the ellipsoidal surface of the 
earth, j^eographic longitude measured east from Greenwich 
in degrees, altitude of the satellite above the oblate 
earth in natuical miles, and geodetic latitude of the 
subvehicle point in degrees. All latitude quantities 
are considered positive north of the equator and nega- 
tive south of the equator. 

20. a, 6, B, A: Right ascension of satellite position, de- 
clination of satellite position, flight path angle, and 
inertial azimuth of velocity vector in units of degrees. 

21. Revolution number, nodal period in minutes, nodal period 
decay rate in minutes per revolution, and nodal regres- 
sion rate in degrees per revolution. Due to the fact 
that the nodal period is calculated by simple subtrac- 
tion of ascending-node crossing times, this quantity 
cannot be determined until two ascending node crossings 
have been detected. Also, since the nodal-period decay 
rate is computed by differencing the nodal-period values 
at successive ascending nodes, this rate cannot be cal- 
culated until three ascending nodes have been crossed. 
The nodal regression rate is computed by differencing 
values of right ascension at successive ascending nodes. 

A . 2 Observation Generation Sample Output 

Figures 20 to 24 contain the output from the observation genera- 
tion test case. There are numbered descriptions corresponding to the 
number codes on the output; however, only the output that has not 
been previously described in Section 4.1 is covered here. Therefore, 
Figures 20 and 21 are included for completeness but need no further 
description (see Section 4.1 for this information). 

1. Bias inputs for data and last END card. 

2. Sigmas. In the case of a data generation run, the 
sigma table is used for the purpose of sepcifying 
the standard deviation of random noise which is to 
be applied to the generated data. Interpretation 
of the printed output shown is otherwise the same 
as described under section for Data Generation* 

3. Station locations, latitude, longitude, and 
altitude. 
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4. Data Limit Specifications. The eight columns from 
left to right contain station identification, data 
interval, minimum elevation angle, maximum eleva- 
tion angle, maximum range in nautical miles, and 
data-generation start and stop times in days, hours, 
and minutes from midnight of epoch day. 

The information shown indicates that the following 
specifications have been given to the program: 

a. Observations are to be generated for station 
AA at 15-second intervals whenever the com- 
puted local elevation angle is above zero 
degrees. 

b. No maximum elevation angle is assumed and no 
maximum range is to be considered (data gen- 
erated for all ranges and for all elevation 
angles greater than zero) . 

c. Station AA is assumed to be active during the 
time interval from epoch until a time 24 hours 
after midnight of epoch day. 

In the case of this particular data generation run, the 
corresponding Specif icat ion- I input items for Stations 
BB through FF are identical to those for Station AA. 
However, this need not be true in general, inasmuch as 
each station is independent of the others with respect 
to these input items . 

5. Parameter indication. The message indicates that the 
range bias parameter has been selected for Station BB. 
In the case of a data generation run, radar parameters 
are selected only for the purpose of applying biases to 
the generated data, in this case a. range bias on Station 
BB data. The message alluding to the correction of para- 
meters therefore should be ignored. 

6. Data Types. These headings define the format of the 
table appearing beneath them. The purpose of this table 
is to define the types of data which are to be generated 
for each station.* Except for the station-identification 
column header, which is printed vertically, the headings 
appear in two horizontal rows, with the top row displaying 
the symbols RANGE on the left and LONG, on the right and 
the second row the symbols SUR.R on the left and U, V on 
the right. The individual headings are interpreted in 
accordance with the following: 
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Heading 
RANGE 
AZMTH 
ELEV. 
R. DOT 
P. DOT 
Q. DOT 
P 

Q 

AZ. DT 
E. DOT 
R. DDT 
MU. VIS 
LAT 
LONG 
SUR. R 

RIGHT 
DOPLR 
LOOK 
VARI 

KAPPA 

ASPCT 
ATT EN 
X, Y, Z 
T-R, D 

HR, ANG 
U, V 



Des c ription or Symbol 

Range 

Local azimuth angle 

Local elevation angle 

Range rate 

Range-rate difference 

Range-rate difference 

Range difference 

Range difference 

Rate of change of local azimuth 

Rate of change of local elevation 

Second time derivative of range 

Mutual visibility 

Latitude of sub-vehicle point 

Longitude of sub-vehicle point 

Surface range, station to sub- 
vehicle point 

Altitude above oblate earth 

Doppler frequency shift 

Look angle 

Variances 



Angle between radius vector 
and local vertical 

Aspect angles 

Signal attenuation 

X, y, £ 

Topocentric right ascension 
and declination 

Topocentric hour angle 

Horizon-scanner angles, u, v 



Units 

n mi 

deg 

deg 

ft/sec 

ft/sec 

ft/sec 

ft 

ft 

deg/sec 

deg/sec 

ft/sec^ 

(indicator) 

de? 

deg 

n mi 

n mi 

cps 

deg 

Same as 

corresp. 

observations 

deg 

deg 
db 

n mi 
deg 

deg 
deg 
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The two horizontal lines associated with each station 
are the line containing the station-identification 
symbol and the line immediately below it. An X on 
the first line of a Y on the second line indicate 
that the corresponding quantity as defined by the table 
header is to be generated. In the present example, the 
indicated data types for Station AA are range, azimuth, 
elevation, range rate, height (altitude) , x, y , z , 
topocentric hour angle, and horizon-sensor angles u 
and V, 

7. Segment transfer messages. 

8. Rise message. The time when the satellite becomes visi- 
ble from a particular station (at the specified minimum 
elevation angle) is obtained by interpolation and printed 
in the manner shown along with the local azimuth angle 
for the correspondinp time. 

9. Time corresponding to generated data. The time shown is 
to be associated with the data quantities appearing on 
the line to the right of the time printout and on the 
line following. The time is given both in hours (0 
through 2A) and minutes of the day as identified at the 
top of the output page, as well as in minutes from start 
(i.e., minutes from eooch) . In the case of the particu- 
lar output shown, time from start corresponds to the 
time of div because the epoch chosen for the run happened 
to be midni'j'bt. 

10. Maximum elevation noint . The time when the elevation 
angle reaches its maximum is obtained by interpolation 
and printed along with the corresponding values for the 
elevation and azimuth angles in the manner shown. 

11. Set message. When the time at which the elevation angle 
reaches the specified minimum from above is obtained by 
interpolation, the print shown at this position occurs. 

12. Duration message. After each pass, a message is orinted 
giving the time in minutes during which the elevation 
angle was above the input minimum value and the range 
was below the input maximum value. 

A . 3 Qj^^t_De termination Sample Output 

This section contains the output from the orbit determination 
test case. Here also, only that output which has not been previously 
described is considered, although all of the output is shown in Fig- 
ures 2 5 through 31. 
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1. The station location card quantities, refraction indi- 
cator, sigma set, latitude, longitude, and altitude. 

2. Each line is the printout of one observation card. 
Given are the date, station ID, both the type and 
set number for these measurements, and the three 
measurements. In this case Line 1 has a Range Rate 
measurement only. 

3. This line contains a Range measurement onlv. 

4. Indicates those parameters chosen in input to be 
estimated. 

5 . In this case , the snherical set of initial condi- 
tions , drag coefficient, GM, J2 and J3 are to be 
estimated. 

6. The number of iterations after which the program 
stops. In this case, one only. 

7. A count of the parameters to be estimated. 

8. A count of the number of observation cards read. 

9. A total obtained by assuming three observations 
per cards. 

10. A count of all the sites entered on station loca- 
tion cards. 

11. The number of locations that all the data has been 
packed into and then stored on scratch tape. 

12. The number of flocks (TF and TR cards) read. 

13. Segment transfer messages. 

14. Differential correction bounds. This sequence 
corresponds to the previously described sequence 
of the parameter indications. 

15. Weighting sigmas. The observation type number 
and the value of the weighting value are shown, 

16. Definition of current orbit. For each iteration 
the current values for the parameters are used 
to compute the Quantities shown 

17. Node print. Each time the integrated trajectory 
crosses the equator (as determined by interpola- 
tion between integration steps) the date, time 
in minutes from midnight of epoch, and the rec- 
tangular elements (x, y, z, x, y, z,) in units 
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of earth radii and earth radii per minute are printed. 
The numerical integration interval for each satellite 
begins at epoch and ends at the time of the latest ob- 
servation for that satellite, A message is printed 
each time any of the events of equator crossing, start 
of thrusting, end of thrusting, or orbit adjust are 
detected during the integration interval. 

18. Station identification and system time of residuals. 
These are to be associated with all residuals which 
appear on the same line. It should be noted that 
system time is defined as seconds from midnight of 
the current day. 

19. Residuals. These are unnormalized differences be- 
tween the input observations (as modified by bias 
or refraction corrections) and corresponding values 
for the same observation types computed from the 
integrated trajectory position at the observation 
times. Up to six residuals for the same time are 
printed on one line. The observation type is indi- 
cated in parentheses iimnediately following the resi- 
dual value in each case. Note that residuals appear 
for the unweighted as well as for the weighted obser- 
vations . 

Identification of the foregoing observation-type in- 
dicators with observation descriptions or symbols 
defined elsewhere in this report is in accordance 
with the following: 

Indicator Description or Symbol Unit 

R Range ft 

A Azimuth deg 

E Elevation deg 

TR Topocentric right ascension deg 

TD Topocentric declination deg 

HA Topocenttic hour angle deg 

GR Geocentric right ascension deg 

GD Geocentric declination deg 

U Horizon scanner in-plane angle deg 

V Horizon scanner cross-plane angle deg 

H Altitude (h) ft 

112. 



Ind i cator Description or Symbol Unit 

X X ft 

Y y ft 

Z z ft 

R Range ft 

P Range difference ft 

C Range difference ft 

RD Range rate ft/sec 

PD Range-rate difference ft/sec 

QD Range-rate difference ft/sec 

20. Observation time. The time identified in Item 2 in terms 
of seconds from midnight is eiven in alternate form* 
wlierein vha dny of Chn mon1:h nnd the hour and minute of 
the dry :.y:c piven r: iuteger:': luv^ jeconclj are gi.veii to 
two cecimr. I »>l;ice , 

21. F>iS summary of residuals. The root-mean-sauare of the 
residuals for each type of observation from each station 
is computed by the residuals editor and the result is 
printed at this location. Included are the station iden- 
tification, number of residuals included in the RMS (i.e., 
the total number for that station and type minus the 
number deleted by the editor on this iteration), RMS, and 
the RMS divided by input weighting sigma. The observa- 
tion-type indicators are the same as those used in the 
residuals print output. Units are feet, degrees and 
seconds. Interpretation of the range-rate summary for 
station MH is that 206 range rate measurements were used 
to compute the rms which equals 1,61216 feet per second. 
RMS over sigma is equal to 26.6209 feet per second. 

22. Iteration number, TRACE-D performs the tracking, or 
orbit determination function, by computing series of 
differential corrections to the parameters selected by 
the user. The iteration number is advanced each time 
the process of computing a set is repeated. The number 
may be interpreted as an indication of the number of 
times the trajectory-integration/ least-squares process 
has been performed. 

23. Observation count. The number of individual observations 
included in the current least-squares computation is 
indicated . 
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24. Convergence indicator. If the weighted RMS for all 
residuals for the current iteration is less than that 
for any previous iteration, then the fitting process 
is converging and the message shown is printed. If 
the RMS obtained on the current iteration is greater 
than the smallest RMS obtained on previous iterations, 
the message 

CURRENT ITERATION IS NOT GOOD 

(RMS = O.xxxxxxxx xx) 

is printed in this position. 

25. Current solution. If the iteration is successful (i.e., 
the overall RMS has been lowered), indicated values are 
the parameter values used in the trajectory, partial 
derivative, and residuals computations which have just 
been completed. In the case of Iteration No. 1, they 
are the input values of the parameters. 

If the iteration is bad, the words 

GO BACK TO 

are printed, and the values of the parameters which so 
far have produced the lowest RMS are recovered from 
memory and printed at this point. Parameter names are 
generally self-explanatory, with the possible exceptions 
of the multiple-vehicle elements. The satellite number 
is printed before the element name for Satellites 2 
through 6. 

26. Current solution in octal digits and machine units. 
These numbers corresponding to those described in Item 
30 are given both in the octal mode and in the units 
used for internal computations. Use of these quantities 
permits bypassing units and number-system conversions 
during input and output. 

27. RMS. This is the quantity which is to be minimized in 
the tracking or orbit determination process and is the 
root-mean-square of the normalized residuals included 
in the least-squares calculations on the current 
iteration. 

28. Corrections. The result of solving the system of normal 
equations associated with the parameter which occupies 
the corresponding position in the current solution block. 
Units are feet, degrees, and seconds. 

29. Bounds. Current values of the numbers used to limit the 
size of corrections. In general, these bounds are 
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automatically increased on a good iteration and auto- 
matically decreased on a bad one. 

30. Bounding indicator. This message will be either 

HITTING BOUNDS 
or 

NOT HITTING BOUNDS 

The first message indicates that the magnitudes of the 
corrections have been controlled by solving the system 
in such a way that the constraint implied by the bounds 
is satisfied, and the latter that the normal equations 
have been solved without applying the bounds. 

31. Next solution. Each value is the sum of the parameter 
value given in the corresponding position under ''current 
solution" and the associated correction. These are the 
parameter values which will be used for the next solution. 

32. Next solution in octal mode and units of earth radii, 
radians, and minutes. 

33. Predicted RMS. If the fitting process is converging in 
a completely linear fashion, this will be the RMS on 
the next iteration. The comparison of this number with 
the current RMS (see Item 32) may be used to measure the 
degree to which the process has already converged. 

34. Sigma of parameters divided by sigma of the normalized 
observations. The numbers given are the square roots 
of the diagonal elements of the inverse normal matrix. 
If certain assumptions are made about the characteris- 
tics of the observation set, the numbers may then be 
taken as the variances on the parameter solutions. 

35. Correlation matrix, correlation coefficients for the 
parameter set. These values are computed directly 
from the covariance matrix (i.e., the inverse normal 
matrix) . Rows and columns are in the same sequence 
as the Item 34 data. 

36. This is the upper* half of the normal matrix. 

37. This is the lower half of the inverse normal matrix. 
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Figure 32. General TRACE-D Flow Chart 
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Figure 33. TRACE-D Segmentation Flow Chart 
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Section 5. 
Programming Information 

5.1 Structures and Flowcharts 

Figures 32 and 33 are both flowcharts of the TRACE-D program; 
however, each is geared to a different level of description. Figure 
32 depicts program paths determined by the function to be performed 
with only basic operations considered, while Figure 33 presents per- 
tinent programming details also. 

Since the TRACE-D program is extremely large (approximately 
22,000 source language cards) and can not be fully contained in core, 
the segmentation mode of operation is employed whereby the program 
becomes several coreloads. Figure 33 shows the levels of segmentation 
employed and the makeup of each and every segment. The fact that only 
the segments along a single vertical path are in core at any one point 
in time solves the storage program. 

The number appearing outside each segment block in Figure 33 
serves to identify its level and position in the structure. Those 
subroutine names underlined are the control programs, which 
among other things, transfer control to the proper segments for each 
major function. 

5.2 Tape Usage 

The TRACE-D program employs many tapes to handle input-output 
data as well as intermediate calculations. The tape requirements for 
each option consists of a basic set common to all functions (listed 
in Table 11) plus a unique option-dependent set. All of the possible 
option-dependent tapes are listed in Tables 12, 13, 14 and those that 
are mandatory are appropriately marked. Please note that no more than 
five tapes from Table 14 may be used in one run (in addition to the 
basic set , of course) . 

The above description has been included for completeness; how- 
ever, as previously mentioned, the user may obtain a running deck 
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for any option he wishes to choose which will set up all appropriate 
trpes for him. 



TAPE NUMBER 



Table 11. 
Basic Set of Tapes 

NAME 
LIBl 
LIB2 



PTAPE* 



FUNCTION 
Program tape //I 
Program tape //2 
System Input Tape 
System Output Tape 
Necessary Scratch Tape 
Necessary Scratch Tape 



* This is a necessary input quantity for all runs. 



TAPE NUMBER 
7 

17 



19 



15 



Table 12. 
Trajectory Prediction Tapes 

NAME FUNCTION 

CTAPE Planetary Coordinate 
Tape 

NTAPE Binary trajectory 

tape (nominal) for 
trajectory differencing 
(program generated) . 

DTAPE Binary difference tape 
for trajectory dif- 
ferencing. 

A binary trajectory 
tape generated by the 
program and saved. 
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TAPE NUMBER 
11 

7 
10 



Table 13. 
Observation Generation Tapes 

NAME 
ETAPE 



FUNCTION 



BCD tape of observations 
program generated. 

CTAPE Planetary tape 

IFLAG(16)* Scratch tape 

(for trajectory) 



* Necessary tape for this option. 







Table lA. 




Orbit 


Determination Tapes 


NUMBER 




NAME 


11 




IBCDI 


12 




IBINI 


7 




CTAPE 


9* 






5* 






10* 







FUNCTION 

BCD tape of observations 

Binary tape of compacted 
observations 

Planetary coordinate 
tape 

Scratch tape (for 
normal matrix) 

Necessary scratch tape 

Scratch tape 
(trajectory) 



* Necessary tapes for this option 
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5 . 3 Subroutines 

This section is devoted to a complete listing of every subroutine: 
in TRACE-D with a functional description and other pertinent facts 
also given. In some cases two or more routines are identical in func- 
tion and differ only slightly in name. This is due to the segmenta- 
tion rule that the same routine name cannot appear in more than one 
logical path (i.e., one complete coreload) . Therefore, since some 
routines were needed in more than one path only the names were changed 
not the content . The routines are listed alphabetically in Table ] S 
and the information given is self-explanatory. 
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5.4 Variable Lists 

In this section are listed several variable arrays used In the 
TRACE-D program. These particular variables were chosen because they 
are the most Important and critical storage areas of TRACE-D. They 
contain all of the Indicators and Input that dictate the logical 
program flow. 
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